Brandvain and Coop. 

Sperm dependent female meiotic 

drive. 



Model I . Taditional female drive 



Set up 

i i I (*Allele and Genotype frequencies*) 
Clear All ["Global^*"] 
f A = 1 - fB; 
f AA = fA"2 + fA fB x; 
f AB = 2 f A fB (1 - x) ; 
f BB = f B " 2 + f A f B x ; 

Drive 

in[7]:= fAA Drive = FullSimplify [fA ( f AA + f AB (1 - d) ) ] ; 

fAB Drive = FullSimplify [f A (f BB + fAB d) + fB (f AA + fAB (1 - d) ) ] ; 
fBB Drive = FullSimplify [fB (fAB d + f BB ) ] ; 

Selection 

in[io]:= wAA = 1; wAB = 1 - hs; wBB = 1 - s; (*genotypic fitnesses*) 

W = FullSimplify [fAA Drive wAA + f AB Drive wAB + fBB Drive wBB] ; (*mean fitness*) 

f AA Sel = FullSimplify [ (f AA Drive * wAA) / w] ; 

f AB Sel = FullSimplify [ (f AB Drive * wAB) / w] ; 

fBB Se i = FullSimplify [ (f BB Drive wBB) / w] ; 

fA Se i = FullSimplify [fAA Se i + fAB Se i / 2] ; 

fB Sel = FullSimplify [fBB Sel + fAB Sel / 2] ; 

Af A = FullSimplify [fA Sel - fA] ; 

AfB = FullSimplify [fB Se i - fB] ; 



Analysis 

For all analyses we assume HWE (i.e. x = 0), we relax this assumption in numerical iterations 
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Invasion 

in[i9]:= scaledDeltaDriveWhenRare = (FullSimplify [Af B / (fB * (1 - fB) ) /. x -> 0] /. fB -> 0) 

1 

out[i9j= - ( - 1 - 2 d ( - 1 + hs ) - hs ) 
2 

Noting that the change in frequency of a standard driver is independent of s when rare, we solve for the 
value of hs required for invasion. 
n[2oj= invasionCrit = Solve [scaledDeltaDriveWhenRare == 0, hs] 

rr -i + 2d 11 

out(20]= -N hs -> \ \ 

11 1 + 2 d J J 

Fixation 

in[2i]:= scaledDeltaDriveWhenCommon = (FullSimplify [AfB / (fB * (1 - fB) ) / . x -» 0] / . fB -> 1) 

l+hs + 2 (d (-1 +hs) - 2 hs + s) 
2 (-1 + S) 

We solve for fixation conditions when recessive or not. 

n[22] = f ixationCritRecessive = Solve [scaledDeltaDriveWhenCommon == 0, s] /. hs -» 0 
f ixationCritNotRecessive = Solve [scaledDeltaDriveWhenCommon == 0, s] 

Out[22]= jjs -> ~ (-1+2 d) JJ 

out[23]= ||s ^ j (-1 + 2 d + 3 hs - 2 d hs) JJ 

Equilibrium 

We identify the equilibrium frequency of the standard driver 

in[24]= eqfB = Solve [ (FullSimplify [ (AfB / . x -> 0) ] ) == 0 , f B] [ [4] ] 

Out[24]= jfB -> 

(8dhs-4ds + A y(-4 (l-2d + hs + 2dhs) (-4hs + 8dhs + 2s-4ds) 
(2 (-4hs + 8dhs + 2s-4ds))J 

We plot this equilibrium frequency assuming a recessive fitness cost. 



- (-8dhs + 4ds) 2 )) j 
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[25]= ContourPlot[If [fB > 1, 1/0, If [fB < 0, 1/0, fB] ] /. eqfB /. hs -> 0, {d, .5, 1} , 

{s, 0, 1}, PlotLegends -» BarLegend [Automatic, LegendLabel -» "fB*"] , PlotLabel -> 

"Equilibrium frequency of a standard driver with a recessive fitness cost", 
FrameLabel -» {"d (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 



Equilibrium frequency of a slandard dii\cr w ith a recessive fitness cos! 




in[26] = FullSimplif y [ f AA Drive + fAB Drive + fBB Drive ] 

Out[26]= 1 

Model 2. Female drive depends on sperm haplotype 
(single pleitropic locus) 

The B allele is transmited with probability, d, in heterozygous females when fertilized by B-bearing 
sperm. 

x represents the deviation from Hardy - Weinberg Equilibrium 



Setup 

m[27]= (*Allele and Genotype frequencies*) 
Clear All ["Global^*"] 
f A = 1 - fB; 
f AA = fA"2 + fA fB x; 
fAB = 2 fA fB (1 - x) ; 
f BB = f B A 2 + f A f B x ; 
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Drive 

in[32]:= (*Genotype frequencies after drive*) 

fAA Drive = FullSimplify [fA ( f AA + fAB / 2) ] ; 

fAB Drive = FullSimplify [fB ( f AA + fAB * (1 -d) ) + f A (fAB / 2 + f BB) ] ; 
fBB Drive = FullSimplify [fB (fAB d + f BB ) ] ; 

Selection 

in[35]:= wAA = 1; wAB = 1 - hs; wBB = 1 - s; (*genotypic fitnesses*) 

W = FullSimplify [fAA Drive wAA + f AB Drive wAB + fBB Drive wBB] ; (*mean fitness*) 

f AA Se i = FullSimplify [ (f AA Drive * wAA) / w] ; 

f AB Sel = FullSimplify [ (f AB Drive * wAB) / w] ; 

fBB Sel = FullSimplify [ (f BB Drive wBB) / w] ; 

fA Sel = FullSimplify [fAA Sel + fAB Sel / 2] ; 

fB Se i = FullSimplify [fBB Se i + fAB Se i / 2] ; 

AfA = FullSimplify [fA Sel - fA] ; 

AfB = FullSimplify [fB Sel - fB] ; 



Analysis 

Note, we assume no deviation from Hardy-Weinberg [i.e. x=0] for all analytical results, and therefore 
these answers are approximations. In the supplamentary material we show thats results of exact recur- 
sions are remarkably consistant from these approximate analystical solutions. 

Assuming the cost of drive is fully recessive [i.e. hs is zero] 
Invasion 

in[44]:= AfBinvade = (FullSimplify [AfB /. hs -> 0 / . x -» 0] / fB"2 /. fB -» 0) 

1 

out[44]= -(-l+d(2-4s)) 
2 

in[45]:= spermDepReceesivelnvade = Solve [AfBinvade == 0, s] 

Ou t[ 45,= {{s^^P}} 

n[46]= plotInvasion4spermDepRecessive = 

Plot[s /. spermDepReceesivelnvade [[1]], {d, .5, 1}, PlotStyle -» {Black, Thick}]; 

[7]= plotRelChange4RarespermDepRecessive = ContourPlot [ {AfBinvade} , {d, 0.5, 1}, 
{s, 0, 1}, PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB"], 
FrameLabel -» { "d (drive)", "s (selection)"}, 

PlotLabel -» "Invasion of recessive sperm dependent driver"]; 
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Fixation 

in[49]:= AfBfix = Ful ISimplif y [ Ful ISimpli f y [ Af B / . hs -» 0 / . x -> 0] / f A] /. f B -» 1 

-1 + 2 d - 2 s 

Out[49]= 

2 - 2 s 

spermDepReceesiveFix = Solve [AfBfix ==0, s] 

Out[50]= {{s -» j (-l + 2d)}} 

in[5i]:= (s /. spermDepReceesiveFix [[1]]) 

1 

outpi]= - ( - 1 + 2 d) 
2 

n[52]= plotFixation4spermDepRecessive = 

Plot[s /. spermDepReceesiveFix [[1]], {d, .5, 1}, PlotStyle -» {Red, Thick}]; 

in[53]= (*Note we artificially rescaled z to be -.1 for all negative values*) 
plotRelChange4CommonSpermDepRecessive = 

ContourPlot[If [s > (s /. spermDepReceesiveFix [ [1] ] ) , - .1, AfBfix] , {d, 0.5, 1} 
{s, 0, 1}, PlotLegends -» BarLegend [Automatic, LegendLabel -» "scaled AB"], 
FrameLabel -» { "d (dirve) " , "s (selection)"}, 

PlotLabel -» "Fixation of recessive sperm dependent driver"]; 
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Bistability Point 

in[55]:= FBbistabSpermDepReceesive = Solve [FullSimplify [AfB / . hs -> 0 / . x -» 0] == 0, fB] [ [4] ] 

r l-2d+4ds, 

Ou.[55]= f fB -» j 

1 -2 s + 4 d s J 

,[56] bistab = ContourPlot[ (If [fB < 0, 0, If [fB > 1, 1, fB] ] ) /. FBbistabSpermDepReceesive, 
{d, .5, 1}, {s, 0, 1}, PlotLegends -» BarLegend [Automatic, LegendLabel -» "fB*"] , 
FrameLabel -» { "d (dirve) " , "s (selection) "}, PlotLabel -» 
"Threshold frequency for fixation of recessive self -promoting driver"]; 
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Assuming the cost of drive is not fully recessive [i.e. hs is 
nonzero] 

Invasion 

Note with any heterozgous cost (i.e. hs > 0) a self - promoting driver cannot invade 
in[58]:= FullSimpli f y [ FullSimpli f y [ Af B / . x -> 0] / f B] / . f B -» 0 

Out[58]= - hS 

Fixation 

AfBfix = FullSimplify [FullSimplify [FullSimplify [AfB / . x -» 0] / f A] / . £B -» 1] 

I + 2 d (-1 + hs) - 3 hs + 2 s 
2 (-1 + S) 

spermDepNotReceesiveFix = Solve [AfBfix == 0 , s] 

||s -> - (-l + 2d + 3hs-2dhs)}} 

spermDepAddFix = Solve [AfBfix == 0 / . hs s / 2 , s] 

II 1 + 2 d J J 



ln[59]:= 
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in[62]:= plotspermDepAddFix = 

Plot[s /. spermDepAddReceesiveFix , {d, .5, 1}, PlotStyle -> {Red, Thick}]; 

ReplaceAII::reps : jspermDepAddReceesiveFixl is neither a list of 

replacement rules nor a valid dispatch table, and so cannot be used for replacing. » 

Bistability Point 

[to]. FBbistabSpermDepNotReceesive = Solve [FullSimplify [AfB /. x -> 0] ==0, fB] [ [3] ] 

Out[63]= jf B -> (-l + 2d + 3hs + 2dhs-4ds- 

^ (-8 hs (-2hs + 4dhs + 2s-4ds) + (l-2d-3hs-2dhs + 4ds) 2 )) / / 
(2 (-2 hs + 4 d hs + 2 s - 4 d s) ) ] 

An Example of a non - recessive driver [Assuming additivity] 

in[64]:= bistab = ContourPlot [ 

(If [fB < 0, 0, If [fB > 1, 1, fB] ] ) /. FBbistabSpermDepNotReceesive /. hs -> (s / 2) , 
{d, .5, 1}, {s, 0, 1}, PlotLegends -» BarLegend [Automatic, LegendLabel -» "fB*"] , 
FrameLabel -» { "d (dirve) " , "s (selection) "}, PlotLabel -» 
"Threshold frequency for fixation of recessive self -promoting driver"]; 

in[65]:= Show[bistab, plotspermDepAddFix] 

Threshold frequency for fixation of recessive self-promoting driver 
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Model 3. Female drive depends on male genotype 
(single pleitropic locus) 

The B allele is transmited with probability, d and dh, in heterozygous females when fertilized BB and AB 
males, respectively. 

x represents the deviation from Hardy - Weinberg Equilibrium 



Setup 

in[68]:= (*Allele and Genotype frequencies*) 
Clear All ["Global^*"] 
f A = 1 - fB; 
f AA = fA"2 + fA fB x; 
f AB = 2 f A fB (1 - x) ; 
f BB = f B A 2 + f A f B x ; 

Drive 

in[73]:= (*Genotype frequencies after drive*) 

fAA Drive = FullSimplif y [ f AA (f AA + fAB / 2) + fAB ( f AA / 2 + f AB (1 - d het ) / 2) ] ; 
fAB Drive = FullSimplify [ 

f AA (fAB / 2 + f BB) + fAB (f AA / 2 + f AB / 2 + f BB (1 - d hom ) ) + f BB (fAA+ fAB / 2) ] ; 
fBB Drive = FullSimplify [fAB (fAB d het / 2 + f BB d hom ) + f BB (fAB / 2 + f BB) ] ; 

Selection 

in[76]:= wAA = 1; wAB = 1 - hs; wBB = 1 - s; (*genotypic fitnesses*) 

W = FullSimplify [fAA Drive wAA + f AB Drive wAB + fBB Drive wBB] ; (*mean fitness*) 

f AA Se i = FullSimplify [ (f AA Drive * wAA) / w] ; 

f AB Sel = FullSimplify [ (f AB Drive * wAB) / w] ; 

fBB Sel = FullSimplify [ (f BB Drive wBB) / w] ; 

fA Sel = FullSimplify [fAA Sel + fAB Sel / 2] ; 

fB Se i = FullSimplify [fBB Se i + fAB Se i / 2] ; 

AfA = FullSimplify [fA Sel - fA] ; 

AfB = FullSimplify [fB Sel - fB] ; 
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in[85]:= FullSimplif y [ f AA Drive + fAB Drive + fBB Drive ] 

Out[85]= 1 



Analysis 

Analytical example - recessive fitness cost 
Invasion 

[86]= invasion4maleDepRecessive = 

Solve[ ( (FullSimplif y[ (AfB /. x -> 0 / . hs -> 0) ] / fB"2) /. fB -> 0) == 0, s] 

° ut[861 ^ {{ s " Vdl?"}' 

in[87]:= plotiInvasion4maleDepRecessive = 

Plot[s /. invasion4maleDepRecessive /. d het -> d hom , 
{dhom, -5, 1}, PlotStyle -> {Black, Thick}]; 

Fixation 

in[88]:= f ixation4maleDepRecessive = 

Solve [ (FullSimplify [ (FullSimplify [ (AfB / . x -> 0 / . hs -> 0) / f A] /. fB -» 1) ]) == 0, s] 



plotFixation4maleDepRecessive = Plot [s /. f ixation4maleDepRecessive /. d het -> d hom 
{dhom, -5, 1}, PlotStyle -> {Red, Thick}]; 

Show[Plot[0, {d hom , 0.5, 1}, AxesLabel -» {"d (drive) assuming d het =d hom =d" , 

"s (selection against drive homozygotes) " } , PlotRange -» { { . 5 , 1}, {0, .5}}, 
PlotLabel -» " Invasion and fixation conditions for a male dependent driver"], 
plotFixation4maleDepRecessive, plotiInvasion4maleDepRecessive] 

Invasion and fixation conditions for a male dependent driver 
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Invasion and fixation conditions for a male dependent drv 




Equlibrium 

in[92]:= FBeqMaleepReceesive = 

Solve [FullSimplif y [AfB / . x -» 0 / . 

r 2 (l-2d hom + 2sd hom ) ] 

Du, [92 ]= f B -> 

L (-1 + 2 s) (-l + 2d hom ) J 



-> 0 / . d het d hoIt 



m[94]= eqfig = ContourPlot [ (If [fB < 0, 0, If [fB > 1, 1, fB] ] ) /. FBeqMaleepReceesive, 

{dhomf 0.5, 1}, {s, 0, .5}, PlotLegends -» BarLegend [Automatic, LegendLabel -> "fB*"] , 
FrameLabel -» { "d (drive) assuming d het =d hom =d" , 

"s (selection, assuming recessive fitness cost)"}]; 
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in[95]:= Show[eqfig, plotiInvasion4maleDepRecessive, plotFixation4maleDepRecessive] 




0.5 0.6 0.7 0.8 0.9 1.0 

d lilnui LisMiining dhet-dhom-d 



Model 4. Female drive depends on sperm haplotype 
(two tightly linked, loci in coupling phase) 

We have one locus with two alleles, A (non-driving) and B (traditional driver), as well as a tightly linked 
locus where one allele modifies drive. Assuming no recombination this functions as a third allele, C. 
This tighly-linked locus is on the B background. 

When C increases drive in heterozgous females it fertilizes, it is a drive enhancer (the B+ allele in our 
ms). 

When C decreases drive in heterozgous females it fertilizes, it is a drive suppressor (the B- allele in our 
ms). 
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Setup 

in[98]:= ClearAll [" Global - *" ] 
fA = . 
f AA = . 
fAB =. 
fAC=. 
fBB =. 
fBC =. 
fCC =. 

minormod = {dl->dO + e} 

(^assuming the sperm acting modifier additively increases drive by epsilon*) ; 
SUMTOONE = { fA -> 1 - (fB + fC) } ; 
HWE = 

{ f AA -> fA"2, fAB -> 2 f A fB, fAC -» 2 f A fC, fBB -» fB"2, fBC -> 2 f B fC, fCC -» fC"2} ; 
GENOFREQS = { f A -» f AA + fAB / 2 + fAC / 2 , 

fB -» fBB + fAB / 2 + fBC / 2 , f C -» f CC + fBC / 2 + fBC / 2} ; 



Drive 

in[ni]:= (*Here we caculate all genotypes after drive. For book- 
keeping purposes we distinguish between reciprocal homozygotes, 
but remove this distinction belowsum them below*) 

ln[112]:= AAn = 

FullSimplif y [ f AA * f AA * 1 + f AA * fAB * 1 / 2 + f AA * fAC * 1 / 2 + f AA * fBB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * (1 - dO) + fAB * fAB * (1 - dO) / 2 + fAB * fAC * (1 - dO) / 2 + 
fAB * fBB * 0 + fAB * fBC * 0 + fAB * f CC * 0 + fAC * f AA * ( 1 - dO) + fAC * fAB * ( 1 - dO) / 2 + 
fAC * fAC * (1 - dO) / 2 + fAC * fBB * 0 + fAC * fBC * 0 + fAC * fCC * 0 + fBB * f AA * 0 + 
fBB * fAB * 0 + fBB * fAC * 0 + fBB * fBB * 0 + fBB * fBC * 0 + fBB * f CC * 0 + fBC * f AA * 0 + 
fBC * fAB * 0 + fBC * fAC * 0 + fBC * fBB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * fBB * o + fee * fBC * o + fee * fee * o + 0] ; 

ABn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 1 / 2 + f AA * fAC * 0 + f AA * fBB * 1 + 
f AA * fBC * 1 / 2 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * ( 1 - dO) / 2 + 
fAB * fAC * 0 + fAB * fBB * (1 - dO) + fAB * fBC * (1 - dO) / 2 + fAB * fCC * 0 + 
fAC * f AA * 0 + fAC * fAB * (1 - dO) / 2 + fAC * fAC * 0 + fAC * fBB * (1 - dO) + 
fAC * fBC * (1 - dO) / 2 + fAC * fCC * 0 + fBB * f AA * 0 + fBB * fAB * 0 + fBB * fAC * 0 + 
fBB * fBB * 0 + fBB * fBC * 0 + fBB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 0 + 
fBC * fAC * 0 + fBC * fBB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * fBB * o + fee * fBC * o + fee * fee * o + 0] ; 

ACn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 1 / 2 + f AA * fBB * 0 + 

f AA * fBC * 1 / 2 + f AA * fCC * 1 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (1 - dl) / 2 + 
fAB * fBB * 0 + fAB * fBC * (1 - dl) / 2 + fAB * fCC * (1 - dl) + fAC * f AA * 0 + fAC * fAB * 0 + 
fAC * fAC * (1 - dl) / 2 + fAC * fBB * 0 + fAC * fBC * (1 - dl) / 2 + f AC * f CC * (1 - dl) + 
fBB * f AA * 0 + fBB * fAB * 0 + fBB * fAC * 0 + fBB * fBB * 0 + fBB * fBC * 0 + fBB * f CC * 0 + 
fBC * f AA * 0 + fBC * fAB * 0 + fBC * fAC * 0 + fBC * fBB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + 

fee * f aa * o + fee * fAB * o + fee * fAC * o + fee * fBB * o + fee * f bc * o + fee * fee * o + 0] ; 

BAn = FullSimplif y [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * fBB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * dO + fAB * fAB * dO / 2 + fAB * fAC * dO / 2 + fAB * fBB * 0 + 
fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + fAC * fBB * 0 + 
fAC * fBC * 0 + fAC * fCC * 0 + fBB * f AA * 1 + fBB * fAB * 1 / 2 + fBB * fAC * 1 / 2 + 
fBB * fBB * 0 + fBB * fBC * 0 + fBB * f CC * 0 + fBC * f AA * 1 / 2 + fBC * fAB * 1 / 4 + 
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fBC * f AC * 1 / 4 + f BC * f BB * 0 + f BC * f BC * 0 + f BC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAc * o + fee * f bb * o + fee * fBc * o + fee * fee * o + 0] ; 

BBn = FullSimplify [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * dO / 2 + fAB * fAC * 0 + fAB * f BB * dO + 
fAB * fBC * dO / 2 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + 
fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 1 / 2 + f BB * fAC * 0 + 
f BB * f BB * 1 + f BB * fBC * 1 / 2 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 1 / 4 + 
fBC * f AC * 0 + f BC * f BB * 1 / 2 + f BC * f BC * 1 / 4 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BCn = FullSimplify [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (dl) / 2 + 

fAB * f BB * 0 + fAB * fBC * (dl) / 2 + fAB * f CC * dl + fAC * f AA * 0 + fAC * fAB * 0 + 

fAC * fAC * 0 + fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + 

f BB * fAC * 1 / 2 + f BB * f BB * 0 + f BB * fBC * 1 / 2 + f BB * f CC * 1 + fBC * f AA * 0 + 

fBC * fAB * 0 + fBC * fAC * 1 / 4 + fBC * f BB * 0 + fBC * fBC * 1 / 4 + fBC * f CC * 1 / 2 + 

fee * f aa * o + fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * f bc * o + fee * fee * o + 0] ; 

CAn = FullSimplify [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + 

fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * dO + fAC * fAB * dO / 2 + fAC * fAC * dO / 2 + 

fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 

f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 1 / 2 + fBC * fAB * 1 / 4 + 

fBC * fAC * 1 / 4 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 1 + 

fCC * fAB * 1 / 2 + fCC * fAC * 1 / 2 + fCC * f BB * 0 + fCC * fBC * 0 + fCC * fCC * 0 + 0] ; 

CBn = FullSimplify [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + fAB * fBC * 0 + 
fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * dO / 2 + fAC * fAC * 0 + fAC * f BB * dO + 
fAC * fBC * dO / 2 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 1 / 4 + 
fBC * fAC * 0 + fBC * f BB * 1 / 2 + fBC * fBC * 1 / 4 + fBC * f CC * 0 + f CC * f AA * 0 + 
fCC * fAB * 1 / 2 + fCC * fAC * 0 + fCC * f BB * 1 + fCC * fBC * 1 / 2 + fCC * fCC * 0 + 0] ; 

CCn = FullSimplify [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + 

fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * dl / 2 + 

fAC * f BB * 0 + fAC * f BC * dl / 2 + fAC * f CC * dl + f BB * f AA * 0 + f BB * f AB * 0 + 

f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 0 + 

fBC * fAC * 1 / 4 + fBC * f BB * 0 + fBC * fBC * 1 / 4 + fBC * f CC * 1 / 2 + f CC * f AA * 0 + 

fCC * fAB * 0 + fCC * fAC * 1 / 2 + fCC * fBB * 0 + fCC * fBC * 1 / 2 + fCC * fCC * 1 + 0] ; 

in[i2i]:= (*Genotype frequencies after drive*) 
fAA Drive = FullSimplify [AAn] ; 
f AB Drive = FullSimplify [ABn + BAn] ; 
f AC Drive = FullSimplify [ACn + CAn] ; 
fBB Drive = FullSimplify [BBn] ; 
fBC Drive = FullSimplify [BCn + CBn] ; 
f CC Drive = FullSimplify [CCn] ; 
(★check, do allele freqs sum to one?*) 
FullSimplify [ 

FullSimplify [f AA Drive + f AB Drive + f AC Drive + fBB Drive + fBC Drive + fCC Drive ] / . HWE / . SUMTOONE] 



Out[127]= 1 
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Selection 



in[i28]= wAA = 1 ; wAC = wAB = 1 - hs; wBB = wBC = wCC = 1 - s; 
W = FullSimplify [ 

(WAA f AA Drive + wAB f AB Drive + wAC f AC Drive + wBB f BB Drive + wBC f BC Drive + wCC f CC Drive ) ] ; 
(♦Because the C allele arises on the B background we assume 
it has the same impact on individual fitness*) 
FullSimplify [w /. HWE /. SUMTOONE /. hs -» 0] 

Out[i30]= 1 + (fB + fC) (-2 dl fC + 2 dO fB (- 1 + fB + fC) - (fB + fC) ( f B + fC - 2 dl f C) ) s 

in[i3i] = f AA Se i = f AA Drive wAA / W; 
f AB Se i = f AB Drive wAB / W; 
f AC Se i = f AC Drive wAC / W; 
f BB Sel = f BB Drive wBB / W; 
f BC Se i = f BC Drive wBC / W; 
f CC Se i = f CC Drive wCC / W; 

fA Se i = FullSimplify [fAA Se i + (fAB Se i + fAC Se i) / 2] ; 
fB Sel = FullSimplify [fBB Sel + (fAB Sel + fBC Se i) / 2] ; 
fC Se i = FullSimplify [fCC Se i + (fAC Se i + fBC Se i) / 2] ; 
Af A = FullSimplify [fA Sel - fA] ; 
AfB = FullSimplify [fB Se i - fB] ; 
AfC = FullSimplify [fCsei - fC] ; 

(*Check: do genotype freqs after selection sum to one?*) 
FullSimplify [fA Se i + fB Se i + fC Se i] 



Note, we assume no deviation from Hardy - Weinberg for all analytical results, and therefore these 
answers are approximations. In the supplamentary material we show thats results of exact recursions 
are remarkably consistant from these approximate analystical solutions. 

Invasion of standard driver - note the driver always invades when it has a recessive fitness cost. 
Reassuringly, this reproduces model 1. 

in[i44] = invasionStandardDriver = Solve [ 

(FullSimplify [ (AfB /. GENOFREQS /. HWE /. SUMTOONE /. fC -> 0) / fB] /. fB -» 0) == 0, hs] 



1 



Analysis 



Analysis - a standard driver [i.e. C is absent] 




-1 + 2 dO 



Fixation of standard driver 
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n[i45]:= f ixationStandardDriver = Solve [ 

(FullSimplify [ (AfB / f A / . GENOFREQS / . HWE / . SUMTOONE / . fC -» 0) ] / . fB -> 1) ==0, s 

y t[ i45]= jjs -> - (-1 + 2 dO + 3 hs - 2 dO hs) J] 

n[i46]:= (*fixation of a standard recessive driver*) 
f ixationStandardDriver / . hs -> 0 

{H; ( - 1 + 2d0) }} 

Equilibrium 

n [i47]= ( *Equilibrium frequency of a standard driver]*) 

n[i48]:= eqfB = Solve [ (AfB / . GENOFREQS / . HWE / . SUMTOONE / . fC -» 0) ==0, fB] [ [4] ] 

ut[i48]= [fB -» 1 8 dO hs - 4 dO s + 

\J -A (1 - 2 dO + hs + 2 dO hs) (-4 hs + 8 dO hs + 2 s - 4 dO s) + (-8 dO hs + 4 dO s) 2 j j 
(2 (-4hs + 8d0hs + 2s-4d0s))} 

n[i49]:= (*Plot of equilibrium frequency of standard driver assuming full recessivity*) 

n[i50]= ContourPlot[If [fB > 1, 1/0, If [fB < 0, 1/0, fB] ] /. eqfB /. hs -> 0, {dO, .5, 1}, 

{s, 0, 1}, PlotLegends -» BarLegend [Automatic, LegendLabel -» "fB*"] , PlotLabel -> 

"Equilibrium frequency of a standard driver with a recessive fitness cost", 
FrameLabel -» {"dO (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 
nquilibi iiim lieqiiencx of a slandurd dn\cv w ith a recessive filness cosl 




Invasion of sperm acting drive modifier tighly linked with the 
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driver, and on the driving background 

(♦change in frequency of the drive modifier when rare and when alleles 
at the drive locus are in drive-viability equilibrium, mutliplied by 
Wbar/fC [this value is always positive and will not influence the sign]*) 
wbarDeltaSpermDrive = FullSimplify [ 

FullSimplify [w AfC / (fC) /. HWE /. SUMTOONE] /. f C -» 0 / . eqfB /. minormod] 
1 

(hs-s) 

2 (1-2 dO) 2 (2 hs-s) 

(-1 - hs - -JY V(l+hs-2d0(2+d0(-2 + s))) (2 hs-s) +2d0 (2 (-1 + dO) (-1 + hs) + s) ) e 



(♦Plotting change this change in frequency when the sperm acting locus is rare, 
increases drive, and when the fitness cost of drive is fully 
recessive. NOTE: This sperm enhancer of drive cannot invade*) 
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in[i53] := ContourPlot [ 

If [fB > 1, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDrive] ] /. eqfB /. e -> 0.01 /. hs 
{d0, 0.5, 1}, {s, 0, 1}, 

PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB + " ] , 
PlotLabel -» "A sperm-acting drive enhancer (the B + haplotype) 

cannot invade\n (assuming a recessive fitness cost of drive)", 
FrameLabel -» {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 

Greater::nord : Invalid comparison with 7.49301 - 83.3298 i attempted. » 
Creater::nord : Invalid comparison with 7.49301 - 83.3298 i attempted. » 
1 

Power::infy : Infinite expression — encountered. » 

0. 

1 

Power::infy : Infinite expression encountered. » 

0. L 

Creater::nord : Invalid comparison with 3.99825 - 83.5704 i attempted. » 
Ceneral::stop : Further output of Creater::nord will be suppressed during this calculation. » 
1 

Power::infy : Infinite expression — encountered. » 

0. 

Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. » 

A sperm-acting drive enhancer (the B : haplotype) cannot invade 
(assuming a recessive fitness cost of time) 




in[i54]:= (*Plotting change this change in frequency when the sperm acting 
locus is rare and decreases drive, and when the fitness cost of 
drive is fully recessive. NOTE: This sperm suppresor always invades*) 
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in[i55]:= ContourPlot [ 

If [fB > 1, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDrive] ] /. eqfB /. e -> -0.01 /. hs -» 0, 
{d0, 0.5, 1}, {s, 0, 1}, 

PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB~ " ] , 
PlotLabel -» "A sperm-acting drive suppressor (the B" haplotype) 

always invades\n (assuming a recessive fitness cost of drive)", 
FrameLabel -» {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 

Creater::nord : Invalid comparison with 7.49301 - 83.3298 i attempted. » 

Creater::nord : Invalid comparison with 7.49301 - 83.3298 i attempted. » 

Power::infy : Infinite express 
Power::infy : Infinite express 

Creater::nord : Invalid comparison with 3.99825 - 83.5704 i attempted. » 
Ceneral::stop : Further output of Creater::nord will be suppressed during this calculation. » 
1 

Power::infy : Infinite expression — encountered. » 

0. 

Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. » 

A sperm-acting drive suppressor (the B~ haplotype) always invades 
(assuming a recessive fitness cost of drive) 
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[156]= (*Plotting change this change in frequency when hte extent 

of drive is mild [dO = .6] and the sperm acting locus is rare 
and increase drive. NOTE: This sperm enhancer never invades*) 
ContourPlot[ (If [fB > 1, 1/0, If [fB < 0.0000001, 0 , wbarDeltaSpermDrive] ] /. eqf B / . 
e -> 0.01 / . d0-> .6), {s, 0, 1}, {hs, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB + " ] , PlotLabel -> 
"A sperm-acting drive enhancer (the B + haplotype) cannot invade\n (assuming 
drive-selection equilibrium of a modest driver, dO = .6)", 
FrameLabel -» { "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 
1 

Power::infy : Infinite expression - encountered. » 

0 

1 

Power::infy : Infinite expression — encountered. » 

0 
1 

Power::infy : Infinite expression — encountered. » 

0 

Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. » 

Creater::nord : Invalid comparison with 3. + 63.2179s attempted. » 

Creater::nord : Invalid comparison with 3. + 63.2179s attempted. » 

A sperm-acting drive enhancer (the B + haplotype) cannot invade 
(assuming drive-seleclion equilibrium of a modes! driver, d() = .6) 



scaled AB + 
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in[i57]:= (*Plotting change this change in frequency when hte extent of 

drive is mild [dO = .98] and the sperm acting locus is rare and 
increase drive. NOTE: This sperm enhancer never invades*) ContourPlot [ 
(If[fB>l, 1/0, If [fB <= 0.001, -0.000000001, wbarDeltaSpermDrive] ] / . eqf B / . 

e -> 0.01 /. dO -> .98) , {s, 0, 1} , {hs, 0, 1} , 
PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB + " ] , PlotLabel -> 
"A sperm-acting drive enhancer (the B + haplotype) cannot invade\n (assuming 
drive-selection equilibrium of a strong driver, dO = .98)", 
FrameLabel -» { "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 
1 

Power::infy : Infinite expression - encountered. » 

0 

Creater::nord : Invalid comparison with 1.02083 - 2.442 i attempted. » 
Creater::nord : Invalid comparison with 1.02083 - 2.442 i attempted. » 
1 

Power::infy : Infinite expression — encountered. » 



Power::infy : Infinite expressio 



tered. ; 



Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. » 

Creater::nord : Invalid comparison with 1.02083 - 1.56808 i attempted. » 

Ceneral::stop : Further output of Creater::nord will be suppressed during this calculation. 

A sperm-acting drive enhancer (the B : haplotype) cannot invade 
(assuming drive-selection equilibrium of a strong driver, dO = .98) 
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Replacement of traditional driver by sperm acting drive 
suppressor tighly linked with the driver, and on the driving 
background 

Equilibrium frequency of drive suppresor 

[158]= eqfC = Solve[ (FullSimplify [Af C /. HWE /. SUMTOONE /. minormod /. fB -> 0]) ==0, fC] [ [4] ] ; 

[159]= (*Equilibrium frequency of a linked, coupled, 

drive suppressor when the fitness costs of drive are recessive*) 

[160]= ContourPlot [If [fC < 1 && fC > 0, (fC) , 1/0] / . eqfC /.hs-»0/.e-» - .001, 

{d0, .5, 1}, {s, 0, 1}, PlotLegends -» BarLegend [Automatic , LegendLabel -> "fB" *"] , 
PlotLabel -» "Equlibrium freq. of sperm-acting drive suppressor (the B" 

haplotype) \n (assuming an initial complete driver, dO = 1)", FrameLabel -» 
{"d0 (drive of the B allele)", "s (selection against drive homozygotes ) " } ] 



tiqulibrium freq. ol sperm-acting drive suppre^or ('.he B haplotype) 
(awiimiiig an initial complete driver, dl) = 1 ) 




[161]= (^Equilibrium frequency of a linked, coupled, 
drive suppressor when drive is complete*) 
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m[i62]:= ContourPlot [ I f [ f C < 1 && fC> 0, (fC) , 1/0] /. eqfC /.hs->0/.e->- .001 /. dO -> 1, 
{s, 0, 1}, {hs, 0, 1}, PlotLegends -> BarLegend [Automatic, LegendLabel -» "fB" *"], 
PlotLabel -> "Equlibrium freq. of a sperm-acting drive suppressor (the 

B" haplotype) \n (assuming an initial complete driver, dO = 1)", 
FrameLabel -» { "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 
1 

Power::infy : Infinite expression — encountered. » 

0 

Less::nord : Invalid comparison with 1.00802 - 2.45019i attempted. » 
Creater::nord : Invalid comparison with 1.00802 - 2.45019 i attempted. » 
Less::nord : Invalid comparison with 1.00802 - 2.45019i attempted. » 
Greater::nord : Invalid comparison with 1.00802 - 2.45019 « attempted. » 
1 

Power::infy : Infinite expression — encountered. » 

0. 
1 

Power::infy : Infinite expression encountered. » 

0. L 

Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. » 
Less::nord : Invalid comparison with 1.00451 - 1.58082 « attempted. » 
Ceneral::stop : Further output of Less::nord will be suppressed during this calculation. » 
Creater::nord : Invalid comparison with 1.00451 - 1.58082 i attempted. » 
Ceneral::stop : Further output of Creater::nord will be suppressed during this calculation. » 
tiqulibrium freq. of a sperm-acting drive suppressor (the 8~ haplotype) 



(assuming an initial complete dri\ er. dO = 1 ) 




Equilibrium frequency of drive/sperm acting suppresor haplotype is often less than that of the 
standard driver it replaces. 
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m[i63]= (*Differecce in equilibrium frequency of the B- and B haplotypes*) 
ContourPlot [ 

If [fC < 1 && fB < 1 && fC> 0 && fB > 0, (fC - fB) , 1 / 0] / . eqfB / . eqfC / . hs -» 0 / . 

e -> - .001, {d0, .5, 1}, {s, 0, 1}, 
PlotLegends -» BarLegend [Automatic, LegendLabel -» "fB" * - fB*"] , PlotLabel -» 
"Difference in equlibrium frequency of a drive suppresor (B" haplotype) \n and 
the driver it replaces (assuming a recessive cost of the drive allele)", 
FrameLabel -» {"d0 (drive of the B allele)", 
"s (selection against drive homozygotes) " } ] 



Difference in equlibrium frequency of a drive suppresor (B haplotype) 
and the driver it replaces (assuming a recessive cosl of the drive allele) 
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n[i64]:= ContourPlot [ 

If [fC < 1 && fB < 1 && fC> 0&& fB > 0, (fC - fB) , 1/0] /. eqfB /. eqfC /. hs -> (s / 4) 

e -> - .001, {d0, .5, 1}, {s, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -» "fB" * - fB*"] , PlotLabel -» 
"Difference in equlibrium frequency of a drive suppresor (B" haplotype) \n 
and the driver it replaces (assuming a heterozygote cost, hs = s/4)", 
FrameLabel -» {"d0 (drive of the B allele)", 
"s (selection against drive homozygotes) " } ] 



Difference in equlibrium frequency of a drive suppresor (B haplotype) 
and the driver it replaces (assuming a hctero/\gole cost, hs = s/4) 




do (drive of the B allele) 

n[i65]= (*change in frequency of a rare traditional driver when the drive- 
sperm suppressor haplotype and the ondriving hplotype 

are at drive selection equilibrium, mutliplied by 
Wbar/fC [this value is always positive and will not influence the sign]*) 
wbarDeltaTradDrive = FullSimplif y [ 

FullSimplify[WAfB/ (fB) /. HWE /. SUMTOONE] /. f B -> 0 / . eqfC /. minormod] ; 

n[i66]:= (*Plotting the change in frequency of a rare traditional 

driving haplotype when the sperm drive suppresor is at drive 
selection balance. ASSUMING the fitness cost of drive is fully 
recessive. NOTE: This sperm enhancer of drive cannot invade*) 
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in[i67] := ContourPlot [ 

If [fC > 1, 1/0, If [fC < 0.0001, 0, wbarDeltaTradDrive] ] / . eqfC / . e -> -0.001 / . hs -» < 
{d0, 0.5, 1}, {s, 0, 1}, 

PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB" ] , 

PlotLabel -> "A fixed drive suppresor (B~ haplotype) cannot be displaced by \n a 

traditional driver (assuming a recessive fitness cost of drive)", 
FrameLabel -» {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 



Power::infy : Infinite expressio 



Creater::nord : Invalid o 
Creater::nord : Invalid o 



with 53.7184 -64.9128 i 
with 53.7184 -64. 9128f 



Power::infy : Infinite expression — eni 
0. 
1 

Power::infy : Infinite expression e 



Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. ; 

A fixed drive suppresor (ZT haplotype) cannot be displaced by 
a traditional drher (assuming a recessive fitness cost of drive) 




in[i68j:= (*Plotting the change in frequency of a rare traditional 

driving haplotype when the sperm drive suppresor is at drive 
selection balance. ASSUMING a fitness cost of drive in homozygotes 
(hs-»s/4) . NOTE: This sperm enhancer of drive cannot invade*) 
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= ContourPlot [ 

If [fC > 1, 1/0, If [fC < 0.0001, 0, wbarDeltaTradDrive]] /. eqfC / . e -> -0.0 

hs -> s / 4, {d0, 0.5, 1}, {s, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AB" ] , 
PlotLabel -» "A fixed sperm-acting drive suppresor cannot be displaced 
by a traditional \n (assuming a recessive fitness cost of drive)", 
FrameLabel -» {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 



Power::infy : Infini 



Creater::nord : Invalid o 
Creater::nord : Invalid o 



with 104.123 - 58.7778« 
with 104.123 - 58.7778f 



Power::infy : Infinite expressio 
Power::infy : Infinite expressio 



encountered. : 



Ceneral::stop : Further output of Power::infy will be suppressed during this calculation. » 
Creater::nord : Invalid comparison with 69.1031 - 97.0611 i attempted. » 
Ceneral::stop : Further output of Creater::nord will be suppressed during this calculation. » 
ContourPlot [ 

If [fC > 1, 1/0, If [fC < 0.0001, 0, wbarDeltaTradDrive]] /. eqfC / . e -» -0.001 /. dO -» 1, 
{s, 0, 1}, {hs, 0, 1}, PlotLegends -» BarLegend [Automatic, LegendLabel -» "scaled AB"] , 
PlotLabel -> "A fixed sperm-acting drive suppresor cannot 

be displaced by a traditional \n (assuming perfect drive)", 
FrameLabel -»{ "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 

A fixed sperm-acting drive suppresor cannot be displaced by a traditional 
(assuming perfect drive) 
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Model 5. Female drive depends on sperm haplotype 
(two tightly linked, loci in repulsion phase) 

We have one locus with two alleles, A (non-driving) and B (traditional driver), as well as a tightly linked 
locus where one allele modifies drive. This tighly-linked locus is on the A background. Assuming no 
recombination this functions as a third allele, C. 

When C increases drive in heterozgous females it fertilizes, it is a drive enhancer (the A+ allele [not 
discussed in our ms]). 

When C decreases drive in heterozgous females it fertilizes, it is a drive suppressor (the A- allele in our 
ms). 



Setup 

Clear All ["Global' *"] 

fA = . 

f AA = . 

fAB =. 

fAC=. 

fBB =. 

fBC =. 

fCC =. 

minormod = {dl-»dO + e} 

(^assuming the sperm acting modifier additively increases drive by epsilon*) ; 
SUMTOONE = { fA -> 1 - (fB + fC) } ; 
HWE = 

{ f AA -> fA"2, fAB -> 2 f A fB, fAC -» 2 f A fC, fBB -» fB~2, fBC -» 2 f B fC, fCC -» fC"2} ; 
GENOFREQS = { f A -» f AA + fAB / 2 + fAC / 2 , 

f B -» fBB + fAB / 2 + fBC / 2 , f C -» f CC + fBC / 2 + fBC / 2 } ; 



Drive 

(*Here we caculate all genotypes after drive. For book- 
keeping purposes we distinguish between reciprocal homozygotes, 
but remove this distinction belowsum them below*) 

AAn = 

FullSimplif y [ f AA * f AA * 1 + f AA * fAB * 1 / 2 + f AA * fAC * 1 / 2 + f AA * fBB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * (1 - dO) + fAB * fAB * (1 - dO) / 2 + fAB * fAC * (1 - dO) / 2 + 
fAB * fBB * 0 + fAB * fBC * 0 + fAB * f CC * 0 + fAC * f AA * 1 / 2 + fAC * fAB * 1 / 4 + 
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f AC * f AC * 1 / 4 + f AC * f BB * 0 + f AC * f BC * 0 + f AC * f CC * 0 + f BB * f AA * 0 + f BB * f AB * 0 + 
f BB * f AC * 0 + f BB * f BB * 0 + f BB * f BC * 0 + f BB * f CC * 0 + f BC * f AA * 0 + f BC * f AB * 0 + 
f BC * f AC * 0 + f BC * f BB * 0 + f BC * f BC * 0 + f BC * f CC * 0 + f CC * f AA * 0 + 
fCC * fAB * 0 + fCC * fAC * 0 + fCC * f BB * 0 + fCC * fBC * 0 + fCC * fCC * 0 + 0] ; 
ABn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 1 / 2 + f AA * fAC * 0 + f AA * f BB * 1 + 

f AA * fBC * 1 / 2 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * ( 1 - dO) / 2 + fAB * fAC * 0 + 

fAB * f BB * (1 - dO) + fAB * fBC * (1 - dO) / 2 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 1 / 4 + 

fAC * f AC * 0 + f AC * f BB * 1 / 2 + f AC * f BC * 1 / 4 + fAC * f CC * 0 + f BB * f AA * 0 + 

f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + 

fBC * fAB * 0 + fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

ACn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 1 / 2 + f AA * f BB * 0 + 

f AA * fBC * 1 / 2 + f AA * fCC * 1 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (1 - dl) / 2 + 
fAB * f BB * 0 + fAB * fBC * (1 - dl) / 2 + fAB * fCC * (1 - dl) + fAC * f AA * 0 + fAC * fAB * 0 + 
fAC * fAC * 1 / 4 + fAC * f BB * 0 + fAC * fBC * 1 / 4 + fAC * f CC * 1 / 2 + f BB * f AA * 0 + 
f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + 
fBC * fAB * 0 + fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BAn = FullSimplif y [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * dO + fAB * fAB * dO / 2 + fAB * fAC * dO / 2 + fAB * f BB * 0 + 
fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + fAC * f BB * 0 + 
fAC * fBC * 0 + fAC * fCC * 0 + f BB * f AA * 1 + f BB * fAB * 1 / 2 + f BB * fAC * 1 / 2 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * dO + fBC * fAB * dO / 2 + 
fBC * fAC * (dO) / 2 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * fCC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BBn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * dO / 2 + fAB * fAC * 0 + fAB * f BB * dO + 
fAB * fBC * dO / 2 + fAB * f CC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + 
fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 1 / 2 + f BB * fAC * 0 + 
f BB * f BB * 1 + f BB * fBC * 1 / 2 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * dO / 2 + 
fBC * fAC * 0 + fBC * f BB * dO + fBC * fBC * (dO) / 2 + fBC * fCC * 0 + fCC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BCn = FullSimplif y [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (dl) / 2 + 

f AB * fBB * 0 + f AB * fBC * (dl) / 2 + f AB * f CC * dl + fAC * f AA * 0 + fAC * fAB * 0 + 

fAC * fAC * 0 + fAC * fBB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + fBB * f AA * 0 + fBB * fAB * 0 + 

fBB * fAC * 1 / 2 + fBB * fBB * 0 + fBB * fBC * 1 / 2 + fBB * f CC * 1 + fBC * f AA * 0 + 

fBC * fAB * 0 + fBC * fAC * (dl) / 2 + fBC * fBB * 0 + fBC * fBC * (dl) / 2 + fBC * fCC * dl + 

fee * f aa * o + fee * fAB * o + fee * fAC * o + fee * fBB * o + fee * f bc * o + fee * fee * o + 0] ; 

CAn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * fBB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * fBB * 0 + 
fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 1 / 2 + fAC * fAB * 1 / 4 + fAC * fAC * 1 / 4 + 
fAC * fBB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + fBB * f AA * 0 + fBB * fAB * 0 + fBB * fAC * 0 + 
fBB * fBB * 0 + fBB * fBC * 0 + fBB * fCC * 0 + fBC * f AA * (1 - dO) + fBC * fAB * (1 - dO) / 2 + 
fBC * fAC * (1 - dO) / 2 + fBC * fBB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + fCC * f AA * 1 + 
fCC * fAB * 1 / 2 + fCC * fAC * 1 / 2 + fCC * fBB * 0 + fCC * fBC * 0 + fCC * fCC * 0 + 0] ; 

CBn = FullSimplif y [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * fBB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * fBB * 0 + fAB * fBC * 0 + 
fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 1 / 4 + fAC * fAC * 0 + fAC * fBB * 1 / 2 + 
fAC * fBC * 1 / 4 + fAC * fCC * 0 + fBB * f AA * 0 + fBB * fAB * 0 + fBB * fAC * 0 + 
fBB * fBB * 0 + fBB * fBC * 0 + fBB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * ( 1 - dO) / 2 + 
fBC * fAC * 0 + fBC * fBB * (1 - dO) + fBC * fBC * (1 - dO) / 2 + fBC * fCC * 0 + f CC * f AA * 0 + 
fCC * fAB * 1 / 2 + fCC * fAC * 0 + fCC * fBB * 1 + fCC * fBC * 1 / 2 + fCC * fCC * 0 + 0] ; 

CCn = FullSimplif y [f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * fBB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * fBB * 0 + fAB * fBC * 0 + 
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f AB * fCC * 0 + f AC * f AA * 0 + f AC * f AB * 0 + f AC * f AC * 1 / 4 + f AC * f BB * 0 + f AC * f BC * 1 / 4 + 

f AC * fCC * 1 / 2 + f BB * f AA * 0 + f BB * f AB * 0 + f BB * f AC * 0 + f BB * f BB * 0 + 

f BB * fBC * 0 + f BB * fCC * 0 + fBC * f AA * 0 + fBC * f AB * 0 + fBC * f AC * (1 - dl) / 2 + 

fBC * f BB * 0 + fBC * fBC * (1 - dl) / 2 + fBC * fCC * (1 - dl) + fCC * f AA * 0 + 

fCC * fAB * 0 + fCC * fAC * 1 / 2 + fCC * f BB * 0 + fCC * fBC * 1 / 2 + fCC * fCC * 1 + 0] ; 

(♦Genotype frequencies after drive*) 

fAA Drive = FullSimplify [AAn] ; 

f AB Drive = FullSimplify [ABn + BAn] ; 

f AC Drive = FullSimplify [ACn + CAn] ; 

fBB Drive = FullSimplify [BBn] ; 

fBC Drive = FullSimplify [BCn + CBn] ; 

fCC Drive = FullSimplify [CCn] ; 

(*check, do allele freqs sum to one?*) 

FullSimplify [ 

FullSimplify [f AA Drive + f AB Drive + f AC Drive + fBB Drive + fBC Drive + fCC Drive ] / . HWE / . SUMTOONE] 

1 



Selection 

wAA = wAC = wCC = 1 ; wAB = wBC = 1 - hs ; wBB = 1 - s ; 
W = FullSimplify [ 

(WAA f AA Drive + wAB f AB Drive + wAC f AC Drive + wBB f BB Drive + wBC f BC Drive + wCC f CC Drive ) ] ; 
(♦Because the C allele arises on the B background we assume 
it has the same impact on individual fitness*) 
FullSimplify [w /. HWE /. SUMTOONE /. hs -» 0] 

1 + (2 dO (-1 + fB) - fB) fB 2 s 

f AA Sel = f AA Drive wAA / W; 
f AB Sel = f AB Drive wAB / W; 
f AC Se i = f AC Drive wAC / W; 
f BB Sel = f BB Drive wBB / W; 
f BCsei = f BC Drive wBC / W; 
fCC S ei = f CC Drive wCC / W; 

fA Sel = FullSimplify [fAA Sel + (fAB Sel + fAC Se i) / 2] ; 
fB Sel = FullSimplify [fBB Sel + (fAB Sel + fBC Se i) / 2] ; 
fC Se i = FullSimplify [fCCsei + (fAC Se i + fBC Se i) / 2] ; 
Af A = FullSimplify [fA Sel - fA] ; 
AfB = FullSimplify [fB Sel - fB] ; 
AfC = FullSimplify [fC Se i - fC] ; 

(*Check: do genotype freqs after selection sum to one?*) 

FullSimplify [fA Sel + fB Sel + fC Se i] 

1 

FullSimplify [wfAC Se i /. HWE /. SUMTOONE] 

2 (-1+ (-1+dO+dl) fB) fC (-1 + fB + fC) 
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Analysis 

Analysis - a standard driver [i.e. C is absent] 

THIS IS THE SAME AS ABOVE, AND WE RENETER IT SIMPLY TO LOAD RESULTS INTO MEM- 
ORY. PLOTS ARE NOT RECREATED 

invasionStandardDriver = Solve [ 

(FullSimplify [ (AfB /. GENOFREQS /. HWE /. SUMTOONE /. fC -» 0) / fB] /. fB -> 0) ==0, 
hs] ; 

f ixationStandardDriver = Solve [ 

(FullSimplify [ (AfB / f A / . GENOFREQS / . HWE / . SUMTOONE / . fC -» 0) ] / . fB -> 1) ==0, 
s] ; 

eqfB = Solve [ (AfB / . GENOFREQS / . HWE / . SUMTOONE / . fC -» 0) ==0, fB] [ [4] ] ; 

Invasion of sperm acting drive modifier tighly linked with the 
driver, on the non-driving background 

(★change in frequency of the drive modifier when rare and when alleles 
at the drive locus are in drive-viability equilibrium, mutliplied by 
Wbar/fC^2 [this value is always positive and will not influence the sign]*) 
wbarDeltaSpermDriveRecessive = FullSimplify [ 

(FullSimplify [WAfC/ (fC"2) /. HWE /. SUMTOONE /. eqfB /. hs -» 0] ) /. minormod] 

(-2 dO s + V2~ V (-1 + 2 dO (2 +d0 (-2 + s) ) ) s ) e 



2 (-1 + 2 dO) s 
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ContourPlot [ 

If [fB > 1, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDriveRecessive] ] /. eqfB /. 

e->-0.01/.hs->0, {d0,0.5, 1}, {s,0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -» "scaled AA~ " ] , 

PlotLabel -> "A sperm-acting drive suppressor linked\n to the nondriving allele 

(the A" haplotype) \n invades when the cost of drive is fully recessive", 
FrameLabel -» {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 



A sperm-acting drive suppressor linked 
to the nondriving allele (the A~ haplots pc i 
invades when the cost of drive is fully recessive 




dO (drive of the B allele) 

(♦change in frequency of the drive modifier when rare and when alleles 
at the drive locus are in drive-viability equilibrium, mutliplied by 
Wbar/fC [this value is always positive and will not influence the sign]*) 

wbarDeltaSpermDrive = FullSimplify [ 

FullSimplify [w AfC / (fC) /. HWE /. SUMTOONE] /. f C -» 0 / . eqfB /. minormod] 

1 hs ( - 1 - hs - 

2 (1 - 2 dO) 2 (2 hs - s) 

VT V(l+hs-2d0 (2+dO (-2 + s))) (2hs-s) +2d0 (2 (-1 + dO) (-1 +hs) + s) ) e 

(★Plotting change this change in frequency when the sperm acting locus is rare, 
increases drive, and when the fitness cost of drive is fully 
recessive. NOTE: This sperm enhancer of drive cannot invade*) 
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ContourPlot |^ 

If [fB > 1, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDrive] ] /. eqfB /. e -> -0.01 /. 

hs -> s / 4, {d0, 0.5, 1}, {s, 0, 1}, 
PlotLegends -» BarLegend [Automatic, LegendLabel -> "scaled AA~ " ] , PlotLabel -» 

"Invasion of a sperm-acting drive suppressor linked to the nondriving\n allele 
(the A" haplotype) , assuming the cost of drive is not fully recessive.", 
FrameLabel -» | "dO (drive of the B allele)", 

"s (selection against driver, assuming het fitness cost, hs= — "Jj 

Invasion of a sperm-acting drive suppressor linked to the nondriving 
allele (the A' haplot\ pe). assuming the com of drive is not iull\ recessive. 




scaled AA~ 

- 0.0004 

- 0.0003 
0.0002 



Model 6. Female drive depends on sperm genotype at 
an unlinked locus 

In addition two our one locus with two alleles, A (non-driving) and B (traditional driver), we have a 
second locus with loci, E and F. We are concerned with the invasion and fixation of the F allele, which 
acts to modify drive when it fertilizes a heterozygote. Drive of the B allele is dO, but changes to d1 when 
fertilized by an 'F bearing sperm 
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Setup 

Clear All ["Global' *"] 

f A = . ; f B = . ; f E = . ; f F = . ; 
f AE = f A f E + D ; 
fAF = fAfF -D; 
f BE = f B f E - D; 
f BF = f B f F + D; 

randomHaps = { f AAEE -» FullSimplif y [ f AE fAE] , 

fABEE -> FullSimplif y [2 fAE f BE] , fBBEE -» FullSimplif y [ f BE f BE] , 

fAAEF -» FullSimplif y [2 fAE fAF] , fABEF -> FullSimplif y [2 (fAE f BF + fAF f BE) ] , 

fBBEF -» FullSimplif y [2 f BE f BF] , fAAFF -» FullSimplif y [ fAF fAF] , 

fABFF -» FullSimplif y [2 fAF f BF] , fBBFF -» FullSimplif y [ fBF f BF] } ; 

sum2one = { f A -» 1 - f B, fE -» 1 - f F} 

minormod = {dl-»dO + e} 
(^assuming the sperm acting modifier additively increases drive by epsilon*) ; 

{ f A -» 1 — f B , f E -» 1 - fF} 

Drive 

fAAEE Drive = FullSimplify [ 

f AAEE * f AAEE * 1 + f AAEE * fABEE * 1 / 2 + f AAEE * fBBEE * 0 + f AAEE * fAAEF * 1 / 2 + 
f AAEE * fABEF * 1 / 4 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + f AAEE * fABFF * 0 + 
f AAEE * fBBFF * 0 + fABEE * f AAEE * (1 - dO) + fABEE * fABEE * (1 - dO) / 2 + fABEE * fBBEE * 0 + 
fABEE * fAAEF * (1 - dO) / 2 + fABEE * fABEF * (1 - dO) / 4 + fABEE * fBBEF * 0 + 
fABEE * fAAFF * 0 + fABEE * fABFF * 0 + fABEE * fBBFF * 0 + fBBEE * f AAEE * 0 + 
fBBEE * fABEE * 0 + fBBEE * fBBEE * 0 + fBBEE * fAAEF * 0 + fBBEE * fABEF * 0 + 
fBBEE * fBBEF * 0 + fBBEE * fAAFF * 0 + fBBEE * fABFF * 0 + fBBEE * fBBFF * 0 + 
fAAEF * f AAEE * 1 / 2 + fAAEF * fABEE * 1 / 4 + fAAEF * fBBEE * 0 + fAAEF * fAAEF * 1 / 4 + 
fAAEF * fABEF * 1 / 8 + fAAEF * fBBEF * 0 + fAAEF * fAAFF * 0 + fAAEF * fABFF * 0 + 



fAAEF 


* fBBFF 


*0 


+ fABEF 


* f AAEE 


* (1 


- dO) / 2 + fABEF * 


fABEE * 


(1 -dO) / 4 + 


fABEF 


* fBBEE 


* 0 


+ fABEF 


* fAAEF 


* (1 


- dO) / 4 + fABEF * 


fABEF * 


(1 -dO) / 8 + 


fABEF 


* fBBEF 


*0 


+ fABEF 


* fAAFF 


*0 + 


fABEF * fABFF * 0 


+ fABEF 


* fBBFF * 0 + 


fBBEF 


* f AAEE 


* 0 


+ fBBEF 


* fABEE 


* 0 + 


fBBEF * fBBEE * 0 


+ fBBEF 


* fAAEF * 0 + 


fBBEF 


* fABEF 


* 0 


+ fBBEF 


* fBBEF 


* 0 + 


fBBEF * fAAFF * 0 


+ fBBEF 


* fABFF * 0 + 


fBBEF 


* fBBFF 


*0 


+ fAAFF 


* f AAEE 


*0 + 


fAAFF * fABEE * 0 


+ fAAFF 


* fBBEE * 0 + 


fAAFF 


* fAAEF 


* 0 


+ fAAFF 


* fABEF 


* 0 + 


fAAFF * fBBEF * 0 


+ fAAFF 


* fAAFF * 0 + 


fAAFF 


* fABFF 


*0 


+ fAAFF 


* fBBFF 


*0 + 


fABFF * f AAEE * 0 


+ fABFF 


* fABEE * 0 + 


fABFF 


* fBBEE 


*0 


+ fABFF 


* fAAEF 


*0 + 


fABFF * fABEF * 0 


+ fABFF 


* fBBEF * 0 + 


fABFF 


* fAAFF 


* 0 


+ fABFF 


* fABFF 


* 0 + 


fABFF * fBBFF * 0 


+ fBBFF 


* f AAEE * 0 + 


fBBFF 


* fABEE 


*0 


+ fBBFF 


* fBBEE 


*0 + 


fBBFF * fAAEF * 0 


+ fBBFF 


* fABEF * 0 + 


fBBFF 


* fBBEF 


* 0 


+ fBBFF 


* fAAFF 


* 0 + 


fBBFF * fABFF * 0 


+ fBBFF 


* fBBFF * 0 + 0] 



f ABEE Drive = FullSimplify [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 1 / 2 + f AAEE * fBBEE * 1 + 
f AAEE * fAAEF * 0 + f AAEE * fABEF * 1 / 4 + f AAEE * fBBEF * 1 / 2 + f AAEE * fAAFF * 0 + 
f AAEE * fABFF * 0 + f AAEE * fBBFF * 0 + fABEE * f AAEE * dO + fABEE * fABEE * 1 / 2 + 
fABEE * fBBEE * (1 - dO) + fABEE * fAAEF * dO / 2 + fABEE * fABEF * 1 / 4 + 
fABEE * fBBEF * (1 - dO) / 2 + fABEE * fAAFF * 0 + fABEE * fABFF * 0 + fABEE * fBBFF * 0 + 
fBBEE * f AAEE * 1 + fBBEE * fABEE * 1 / 2 + fBBEE * fBBEE * 0 + fBBEE * fAAEF * 1 / 2 + 
fBBEE * fABEF * 1 / 4 + fBBEE * fBBEF * 0 + fBBEE * fAAFF * 0 + fBBEE * fABFF * 0 + 
fBBEE * fBBFF * 0 + fAAEF * f AAEE * 0 + fAAEF * fABEE * 1 / 4 + fAAEF * fBBEE * 1 / 2 + 
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f AAEF * f AAEF * 0 + f AAEF * f ABEF * 1 / 8 + f AAEF * f BBEF * 1 / 4 + f AAEF * f AAFF * 0 + 

f AAEF * f ABFF * 0 + f AAEF * f BBFF * 0 + f ABEF * f AAEE * dO / 2 + f ABEF * f ABEE * 1 / 4 + 

f ABEF * fBBEE * (1 - dO) / 2 + f ABEF * f AAEF * dO / 4 + f ABEF * f ABEF * 1 / 8 + 

f ABEF * f BBEF * (1 - dO) / 4 + f ABEF * f AAFF * 0 + f ABEF * f ABFF * 0 + f ABEF * f BBFF * 0 + 

f BBEF * f AAEE * 1 / 2 + f BBEF * f ABEE * 1 / 4 + f BBEF * fBBEE * 0 + f BBEF * f AAEF * 1 / 4 + 

f BBEF * f ABEF * 1 / 8 + f BBEF * f BBEF * 0 + f BBEF * f AAFF * 0 + f BBEF * f ABFF * 0 + 

f BBEF * f BBFF * 0 + f AAFF * f AAEE * 0 + f AAFF * f ABEE * 0 + f AAFF * fBBEE * 0 + 

f AAFF * f AAEF * 0 + f AAFF * f ABEF * 0 + f AAFF * f BBEF * 0 + f AAFF * f AAFF * 0 + 

f AAFF * f ABFF * 0 + f AAFF * f BBFF * 0 + f ABFF * f AAEE * 0 + f ABFF * f ABEE * 0 + 

f ABFF * fBBEE * 0 + f ABFF * f AAEF * 0 + f ABFF * f ABEF * 0 + f ABFF * f BBEF * 0 + 

f ABFF * f AAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 

f BBFF * f ABEE * 0 + f BBFF * fBBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * f ABEF * 0 + 

f BBFF * f BBEF * 0 + f BBFF * f AAFF * 0 + f BBFF * f ABFF * 0 + f BBFF * f BBFF * 0 + 0] 

f BBEE Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * f ABEE * 0 + f AAEE * fBBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 0 + f AAEE * f BBEF * 0 + f AAEE * f AAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * f BBFF * 0 + f ABEE * f AAEE * 0 + f ABEE * f ABEE * dO / 2 + 
f ABEE * fBBEE * dO + fABEE * f AAEF * 0 + fABEE * f ABEF * (dO / 4) + 
f ABEE * f BBEF * dO / 2 + fABEE * f AAFF * 0 + fABEE * f ABFF * 0 + fABEE * f BBFF * 0 + 
fBBEE * f AAEE * 0 + fBBEE * fABEE * 1 / 2 + fBBEE * fBBEE * 1 + fBBEE * f AAEF * 0 + 
fBBEE * f ABEF * 1 / 4 + fBBEE * f BBEF * 1 / 2 + fBBEE * f AAFF * 0 + fBBEE * f ABFF * 0 + 
fBBEE * f BBFF * 0 + f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * fBBEE * 0 + 
f AAEF * f AAEF * 0 + f AAEF * f ABEF * 0 + f AAEF * f BBEF * 0 + f AAEF * f AAFF * 0 + 
f AAEF * f ABFF * 0 + f AAEF * f BBFF * 0 + f ABEF * f AAEE * 0 + f ABEF * fABEE * (dO / 4) + 
f ABEF * fBBEE * dO / 2 + f ABEF * f AAEF * 0 + f ABEF * f ABEF * (dO / 8) + 
f ABEF * f BBEF * dO / 4 + f ABEF * f AAFF * 0 + f ABEF * f ABFF * 0 + f ABEF * f BBFF * 0 + 
f BBEF * f AAEE * 0 + f BBEF * fABEE * 1 / 4 + f BBEF * fBBEE * 1 / 2 + f BBEF * f AAEF * 0 + 
f BBEF * f ABEF * 1 / 8 + f BBEF * f BBEF * 1 / 4 + f BBEF * f AAFF * 0 + f BBEF * f ABFF * 0 + 
f BBEF * f BBFF * 0 + f AAFF * f AAEE * 0 + f AAFF * fABEE * 0 + f AAFF * fBBEE * 0 + 
f AAFF * f AAEF * 0 + f AAFF * f ABEF * 0 + f AAFF * f BBEF * 0 + f AAFF * f AAFF * 0 + 
f AAFF * f ABFF * 0 + f AAFF * f BBFF * 0 + f ABFF * f AAEE * 0 + f ABFF * fABEE * 0 + 
f ABFF * fBBEE * 0 + f ABFF * f AAEF * 0 + f ABFF * f ABEF * 0 + f ABFF * f BBEF * 0 + 
f ABFF * f AAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 
f BBFF * fABEE * 0 + f BBFF * fBBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * f ABEF * 0 + 
f BBFF * f BBEF * 0 + f BBFF * f AAFF * 0 + f BBFF * f ABFF * 0 + f BBFF * f BBFF * 0 + 0] 

f AAEF Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * fBBEE * 0 + 
f AAEE * f AAEF * 1 / 2 + f AAEE * f ABEF * 1 / 4 + f AAEE * f BBEF * 0 + f AAEE * f AAFF * 1 + 
f AAEE * f ABFF * 1 / 2 + f AAEE * f BBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * fBBEE * 0 + fABEE * f AAEF * (1 - dl) / 2 + fABEE * f ABEF * (1 - dl) / 4 + 
fABEE * f BBEF * 0 + fABEE * f AAFF * (1 - dl) + fABEE * f ABFF * (1 - dl) / 2 + 
fABEE * f BBFF * 0 + fBBEE * f AAEE * 0 + fBBEE * fABEE * 0 + fBBEE * fBBEE * 0 + 
fBBEE * f AAEF * 0 + fBBEE * f ABEF * 0 + fBBEE * f BBEF * 0 + fBBEE * f AAFF * 0 + 
fBBEE * f ABFF * 0 + fBBEE * f BBFF * 0 + f AAEF * f AAEE * 1 / 2 + f AAEF * fABEE * 1 / 4 + 
f AAEF * fBBEE * 0 + f AAEF * f AAEF * 1 / 2 + f AAEF * f ABEF * 1 / 4 + f AAEF * f BBEF * 0 + 
f AAEF * f AAFF * 1 / 2 + f AAEF * f ABFF * 1 / 4 + f AAEF * f BBFF * 0 + f ABEF * f AAEE * (1 - dO) / 2 
f ABEF * fABEE * (1 - dO) / 4 + f ABEF * fBBEE * 0 + f ABEF * f AAEF * ( (1 - dO) / 4 + (1 - dl) / 4) 
f ABEF * f ABEF * ( (1 - dO) / 8 + (1 - dl) / 8) + f ABEF * f BBEF * 0 + f ABEF * f AAFF * (1 - dl) / 2 
f ABEF * f ABFF * (1 - dl) / 4 + f ABEF * f BBFF * 0 + f BBEF * f AAEE * 0 + f BBEF * fABEE * 0 + 
f BBEF * fBBEE * 0 + f BBEF * f AAEF * 0 + f BBEF * f ABEF * 0 + f BBEF * f BBEF * 0 + 
f BBEF * f AAFF * 0 + f BBEF * f ABFF * 0 + f BBEF * f BBFF * 0 + f AAFF * f AAEE * 1 + 
f AAFF * fABEE * 1 / 2 + f AAFF * fBBEE * 0 + f AAFF * f AAEF * 1 / 2 + f AAFF * f ABEF * 1 / 4 + 
f AAFF * f BBEF * 0 + f AAFF * f AAFF * 0 + f AAFF * f ABFF * 0 + f AAFF * f BBFF * 0 + 
f ABFF * f AAEE * dO + f ABFF * fABEE * (1 - dO) / 2 + f ABFF * fBBEE * 0 + 
f ABFF * f AAEF * (1 - dO) / 2 + f ABFF * f ABEF * (1 - dO) / 4 + f ABFF * f BBEF * 0 + 
f ABFF * f AAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 
f BBFF * fABEE * 0 + f BBFF * fBBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * f ABEF * 0 + 
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fBBFF * fBBEF * 0 + fBBFF * fAAFF * 0 + fBBFF * f ABFF * 0 + fBBFF * fBBFF * 0 + 0] 
f ABEF Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 1 / 4 + f AAEE * fBBEF * 1 / 2 + f AAEE * fAAFF * 0 + 
f AAEE * f ABFF * 1 / 2 + f AAEE * fBBFF * 1 + f ABEE * f AAEE * 0 + f ABEE * f ABEE * 0 + 
f ABEE * f BBEE * 0 + f ABEE * f AAEF * dl / 2 + f ABEE * f ABEF * 1 / 4 + 
fABEE * fBBEF * (1 - dl) / 2 + fABEE * fAAFF * dl + fABEE * f ABFF * 1 / 2 + 
f ABEE * fBBFF * (1 - dl) + f BBEE * f AAEE * 0 + f BBEE * fABEE * 0 + f BBEE * f BBEE * 0 + 
f BBEE * f AAEF * 1 / 2 + f BBEE * f ABEF * 1 / 4 + f BBEE * fBBEF * 0 + f BBEE * fAAFF * 1 + 
f BBEE * f ABFF * 1 / 2 + f BBEE * fBBFF * 0 + f AAEF * f AAEE * 0 + f AAEF * fABEE * 1 / 4 + 
f AAEF * f BBEE * 1 / 2 + f AAEF * f AAEF * 0 + f AAEF * f ABEF * 1 / 4 + f AAEF * fBBEF * 1 / 2 + 
f AAEF * fAAFF * 0 + f AAEF * f ABFF * 1 / 4 + f AAEF * fBBFF * 1 / 2 + f ABEF * f AAEE * dO / 2 + 
f ABEF * fABEE * 1 / 4 + f ABEF * f BBEE * (1 - dO) / 2 + f ABEF * f AAEF * (dO / 4 + dl / 4) + 
f ABEF * f ABEF * 1 / 4 + f ABEF * fBBEF * ( (1 - dO) / 4 + (1 - dl) / 4) + f ABEF * fAAFF * dl / 2 + 
f ABEF * f ABFF * 1 / 4 + f ABEF * fBBFF * (1 - dl) / 2 + fBBEF * f AAEE * 1 / 2 + 
fBBEF * fABEE * 1 / 4 + fBBEF * f BBEE * 0 + fBBEF * f AAEF * 1 / 2 + fBBEF * f ABEF * 1 / 4 + 
fBBEF * fBBEF * 0 + fBBEF * fAAFF * 1 / 2 + fBBEF * f ABFF * 1 / 4 + fBBEF * fBBFF * 0 + 
fAAFF * f AAEE * 0 + fAAFF * fABEE * 1 / 2 + fAAFF * f BBEE * 1 + fAAFF * f AAEF * 0 + 
fAAFF * f ABEF * 1 / 4 + fAAFF * fBBEF * 1 / 2 + fAAFF * fAAFF * 0 + fAAFF * f ABFF * 0 + 
fAAFF * fBBFF * 0 + f ABFF * f AAEE * (1 - dO) + f ABFF * fABEE * 1 / 2 + f ABFF * f BBEE * dO + 
f ABFF * f AAEF * dO / 2 + f ABFF * f ABEF * 1 / 4 + f ABFF * fBBEF * (1 - dO) / 2 + 
f ABFF * fAAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * fBBFF * 0 + fBBFF * f AAEE * 1 + 
fBBFF * fABEE * 1 / 2 + fBBFF * f BBEE * 0 + fBBFF * f AAEF * 1 / 2 + fBBFF * f ABEF * 1 / 4 + 
fBBFF * fBBEF * 0 + fBBFF * fAAFF * 0 + fBBFF * f ABFF * 0 + fBBFF * fBBFF * 0 + 0] 
f BBEF Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 0 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * fBBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * f BBEE * 0 + fABEE * f AAEF * 0 + fABEE * f ABEF * dl / 4 + fABEE * fBBEF * dl / 2 + 
fABEE * fAAFF * 0 + fABEE * f ABFF * dl / 2 + fABEE * fBBFF * dl + f BBEE * f AAEE * 0 + 
f BBEE * fABEE * 0 + f BBEE * f BBEE * 0 + f BBEE * f AAEF * 0 + f BBEE * f ABEF * 1 / 4 + 
f BBEE * fBBEF * 1 / 2 + f BBEE * fAAFF * 0 + f BBEE * f ABFF * 1 / 2 + f BBEE * fBBFF * 1 + 
f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * f BBEE * 0 + f AAEF * f AAEF * 0 + 
f AAEF * f ABEF * 0 + f AAEF * fBBEF * 0 + f AAEF * fAAFF * 0 + f AAEF * f ABFF * 0 + 
f AAEF * fBBFF * 0 + f ABEF * f AAEE * 0 + f ABEF * fABEE * (dO / 4) + f ABEF * f BBEE * (dO / 2) + 
f ABEF * f AAEF * 0 + f ABEF * f ABEF * (dO / 8 + dl / 8) + f ABEF * fBBEF * (dO / 4 + dl / 4) + 
f ABEF * fAAFF * 0 + f ABEF * f ABFF * dl / 4 + f ABEF * fBBFF * (dl / 2) + 
fBBEF * f AAEE * 0 + fBBEF * fABEE * 1 / 4 + fBBEF * f BBEE * 1 / 2 + fBBEF * f AAEF * 0 + 
fBBEF * f ABEF * 1 / 4 + fBBEF * fBBEF * 1 / 2 + fBBEF * fAAFF * 0 + fBBEF * f ABFF * 1 / 4 + 
fBBEF * fBBFF * 1 / 2 + fAAFF * f AAEE * 0 + fAAFF * fABEE * 0 + fAAFF * f BBEE * 0 + 
fAAFF * f AAEF * 0 + fAAFF * f ABEF * 0 + fAAFF * fBBEF * 0 + fAAFF * fAAFF * 0 + 
fAAFF * f ABFF * 0 + fAAFF * fBBFF * 0 + f ABFF * f AAEE * 0 + f ABFF * fABEE * dO / 2 + 
f ABFF * f BBEE * (1 - dO) + f ABFF * f AAEF * 0 + f ABFF * f ABEF * dO / 4 + f ABFF * fBBEF * dO / 2 
f ABFF * fAAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * fBBFF * 0 + fBBFF * f AAEE * 0 + 
fBBFF * fABEE * 1 / 2 + fBBFF * f BBEE * 1 + fBBFF * f AAEF * 0 + fBBFF * f ABEF * 1 / 4 + 
fBBFF * fBBEF * 1 / 2 + fBBFF * fAAFF * 0 + fBBFF * f ABFF * 0 + fBBFF * fBBFF * 0 + 0] 
f AAFF Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 0 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * fBBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * f BBEE * 0 + fABEE * f AAEF * 0 + fABEE * f ABEF * 0 + fABEE * fBBEF * 0 + 
fABEE * fAAFF * 0 + fABEE * f ABFF * 0 + fABEE * fBBFF * 0 + f BBEE * f AAEE * 0 + 
f BBEE * fABEE * 0 + f BBEE * f BBEE * 0 + f BBEE * f AAEF * 0 + f BBEE * f ABEF * 0 + 
f BBEE * fBBEF * 0 + f BBEE * fAAFF * 0 + f BBEE * f ABFF * 0 + f BBEE * fBBFF * 0 + 
f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * f BBEE * 0 + f AAEF * f AAEF * 1 / 4 + 
f AAEF * f ABEF * 1 / 8 + f AAEF * fBBEF * 0 + f AAEF * fAAFF * 1 / 2 + 
f AAEF * f ABFF * 1 / 4 + f AAEF * fBBFF * 0 + f ABEF * f AAEE * 0 + f ABEF * fABEE * 0 + 
f ABEF * f BBEE * 0 + f ABEF * f AAEF * (1 - dl) / 4 + f ABEF * f ABEF * (1 - dl) / 8 + 
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fABEF * fBBEF * 0 + fABEF * fAAFF * (1 - dl) / 2 + f ABEF * fABFF * (1 - dl) / 4 + 

f ABEF * f BBFF * 0 + fBBEF * f AAEE * 0 + fBBEF * f ABEE * 0 + fBBEF * f BBEE * 0 + 

fBBEF * f AAEF * 0 + fBBEF * fABEF * 0 + fBBEF * fBBEF * 0 + fBBEF * fAAFF * 0 + 

fBBEF * fABFF * 0 + fBBEF * f BBFF * 0 + fAAFF * f AAEE * 0 + fAAFF * f ABEE * 0 + 

fAAFF * f BBEE * 0 + fAAFF * f AAEF * 1 / 2 + fAAFF * fABEF * 1 / 4 + fAAFF * fBBEF * 0 + 

fAAFF * fAAFF * 1 + fAAFF * fABFF * 1 / 2 + fAAFF * f BBFF * 0 + fABFF * f AAEE * 0 + 

fABFF * fABEE * 0 + fABFF * f BBEE * 0 + fABFF * f AAEF * (1 - dl) / 2 + 

fABFF * fABEF * (1 - dl) / 4 + fABFF * fBBEF * 0 + fABFF * fAAFF * (1 - dl) + 

fABFF * fABFF * (1 - dl) / 2 + fABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 

f BBFF * fABEE * 0 + f BBFF * f BBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * fABEF * 0 + 

f BBFF * fBBEF * 0 + f BBFF * fAAFF * 0 + f BBFF * fABFF * 0 + f BBFF * f BBFF * 0 + 0] 

f ABFF Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * fABEF * 0 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + 
f AAEE * fABFF * 0 + f AAEE * f BBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * f BBEE * 0 + fABEE * f AAEF * 0 + fABEE * fABEF * 0 + fABEE * fBBEF * 0 + 
fABEE * fAAFF * 0 + fABEE * fABFF * 0 + fABEE * f BBFF * 0 + f BBEE * f AAEE * 0 + 
f BBEE * fABEE * 0 + f BBEE * f BBEE * 0 + f BBEE * f AAEF * 0 + f BBEE * fABEF * 0 + 
f BBEE * fBBEF * 0 + f BBEE * fAAFF * 0 + f BBEE * fABFF * 0 + f BBEE * f BBFF * 0 + 
f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * f BBEE * 0 + f AAEF * f AAEF * 0 + 
f AAEF * fABEF * 1 / 8 + f AAEF * fBBEF * 1 / 4 + f AAEF * fAAFF * 0 + f AAEF * fABFF * 1 / 4 + 
f AAEF * f BBFF * 1 / 2 + fABEF * f AAEE * 0 + fABEF * fABEE * 0 + fABEF * f BBEE * 0 + 
fABEF * f AAEF * dl / 4 + fABEF * fABEF * 1 / 8 + fABEF * fBBEF * (1 - dl) / 4 + 
fABEF * fAAFF * dl / 2 + fABEF * fABFF * 1 / 4 + fABEF * f BBFF * (1 - dl) / 2 + 
fBBEF * f AAEE * 0 + fBBEF * fABEE * 0 + fBBEF * f BBEE * 0 + fBBEF * f AAEF * 1 / 4 + 
fBBEF * fABEF * 1 / 8 + fBBEF * fBBEF * 0 + fBBEF * fAAFF * 1 / 2 + fBBEF * fABFF * 1 / 4 + 
fBBEF * f BBFF * 0 + fAAFF * f AAEE * 0 + fAAFF * fABEE * 0 + fAAFF * f BBEE * 0 + 
fAAFF * f AAEF * 0 + fAAFF * fABEF * 1 / 4 + fAAFF * fBBEF * 1 / 2 + fAAFF * fAAFF * 0 + 
fAAFF * fABFF * 1 / 2 + fAAFF * f BBFF * 1 + fABFF * f AAEE * 0 + fABFF * fABEE * 0 + 
fABFF * f BBEE * 0 + fABFF * f AAEF * dl / 2 + fABFF * fABEF * 1 / 4 + fABFF * fBBEF * (1 - dl) / 2 + 
fABFF * fAAFF * dl + fABFF * fABFF * 1 / 2 + fABFF * f BBFF * (1 - dl) + f BBFF * f AAEE * 0 + 
f BBFF * fABEE * 0 + f BBFF * f BBEE * 0 + f BBFF * f AAEF * 1 / 2 + f BBFF * fABEF * 1 / 4 + 
f BBFF * fBBEF * 0 + f BBFF * fAAFF * 1 + f BBFF * fABFF * 1 / 2 + f BBFF * f BBFF * 0 + 0] 

f BBFF Drive = FullSimplif y [ f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * fABEF * 0 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + 
f AAEE * fABFF * 0 + f AAEE * f BBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * f BBEE * 0 + fABEE * f AAEF * 0 + fABEE * fABEF * 0 + fABEE * fBBEF * 0 + 
fABEE * fAAFF * 0 + fABEE * fABFF * 0 + fABEE * f BBFF * 0 + f BBEE * f AAEE * 0 + 
f BBEE * fABEE * 0 + f BBEE * f BBEE * 0 + f BBEE * f AAEF * 0 + f BBEE * fABEF * 0 + 
f BBEE * fBBEF * 0 + f BBEE * fAAFF * 0 + f BBEE * fABFF * 0 + f BBEE * f BBFF * 0 + 
f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * f BBEE * 0 + f AAEF * f AAEF * 0 + 
f AAEF * fABEF * 0 + f AAEF * fBBEF * 0 + f AAEF * fAAFF * 0 + f AAEF * fABFF * 0 + 
f AAEF * f BBFF * 0 + fABEF * f AAEE * 0 + fABEF * fABEE * 0 + fABEF * f BBEE * 0 + 
fABEF * f AAEF * 0 + fABEF * fABEF * dl / 8 + fABEF * fBBEF * dl / 4 + 
fABEF * fAAFF * 0 + fABEF * fABFF * dl / 4 + fABEF * f BBFF * dl / 2 + 
fBBEF * f AAEE * 0 + fBBEF * fABEE * 0 + fBBEF * f BBEE * 0 + fBBEF * f AAEF * 0 + 
fBBEF * fABEF * 1 / 8 + fBBEF * fBBEF * 1 / 4 + fBBEF * fAAFF * 0 + fBBEF * fABFF * 1 / 4 + 
fBBEF * f BBFF * 1 / 2 + fAAFF * f AAEE * 0 + fAAFF * fABEE * 0 + fAAFF * f BBEE * 0 + 
fAAFF * f AAEF * 0 + fAAFF * fABEF * 0 + fAAFF * fBBEF * 0 + fAAFF * fAAFF * 0 + 
fAAFF * fABFF * 0 + fAAFF * f BBFF * 0 + fABFF * f AAEE * 0 + fABFF * fABEE * 0 + 
fABFF * f BBEE * 0 + fABFF * f AAEF * 0 + fABFF * fABEF * dl / 4 + fABFF * fBBEF * dl / 2 + 
fABFF * fAAFF * 0 + fABFF * fABFF * dl / 2 + fABFF * f BBFF * (dl) + f BBFF * f AAEE * 0 + 
f BBFF * fABEE * 0 + f BBFF * f BBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * fABEF * 1 / 4 + 
f BBFF * fBBEF * 1 / 2 + f BBFF * fAAFF * 0 + f BBFF * fABFF * 1 / 2 + f BBFF * f BBFF * 1 + 0] 
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- (4 fAAEE + 2 (fAAEF + fABEE) + fABEF) ( 2 f AAEE + f AAEF - ( - 1 + dO ) (2 fABEE + fABEF) ) 
8 

1 

- ( fAAEF (2 fABEE + 4 dO fABEE + fABEF + 2 dO fABEF + 8 f BBEE + 4 f BBEF) + 
8 

2 fAAEE ( (2 + 4 dO) fABEE + fABEF + 2 dO fABEF + 8 f BBEE + 4 f BBEF ) + 
(2 fABEE + fABEF) (2 fABEE + fABEF - (-3 + 2 dO) (2 f BBEE + f BBEF ) ) ) 

1 

- (dO (2 fABEE + fABEF) + 2 f BBEE + f BBEF) (2 fABEE + fABEF + 4 f BBEE + 2 f BBEF) 



- ( 4 fAAEF 2 + 2 fAAEE (4 fAAEF + 8 fAAFF + 3 fABEF - 2 dO fABEF + 2 fABFF + 4 dO fABFF) - 

8 V 

(2 fABEE + fABEF) ((-6 + 4dl) f AAFF + (-2+dO+dl) ( fABEF + 2 fABFF )) + 

fAAEF (8 fAAFF - 2 ( (-3 + 2 dl) fABEE + (-3 + dO + dl) fABEF + (-3 + 2 dO) fABFF) ) ) 

1 

- (2 fABEE fABEF + fABEF 2 + 4 fABEE fABFF + 2 fABEF fABFF + 3 fABEF f BBEE - 
4 

2 dO fABEF f BBEE + 2 fABFF f BBEE + 4 dO fABFF f BBEE + 3 fABEE f BBEF - 
2 dl fABEE f BBEF + 3 fABEF f BBEF - dO fABEF f BBEF - dl fABEF f BBEF + 3 fABFF f BBEF - 
2 dO fABFF f BBEF + fAAFF ( (2 + 4 dl) fABEE + fABEF + 2 dl fABEF + 8 f BBEE + 4 f BBEF) - 
( - 3 + 2 dl ) (2 fABEE + fABEF ) f BBFF + 

fAAEE ( fABEF + 2 dO fABEF + 6 fABFF - 4 dO fABFF + 4 f BBEF + 8 f BBFF ) + fAAEF 
(fABEE + 2 dl fABEE + (1 + dO + dl) fABEF + fABFF + 2 dO fABFF + 4 ( f BBEE + f BBEF + f BBFF) ) ) 



- (dO (fABEF 2 + 2 fABEE (fABEF + 2 fABFF) + 

8 X X 

4 fABFF (-2 f BBEE + f BBEF) + 2 fABEF (fABFF + 2 f BBEE + f BBEF ) ) + 
2 (fABFF (6 f BBEE + f BBEF) + fABEF ( f BBEE + f BBEF + f BBFF ) + 

(fABEE + 4 f BBEE + 2 f BBEF) ( f BBEF + 2 f BBFF) ) + 
dl (2 fABEE + fABEF) (fABEF + 2 (fABFF + f BBEF + 2 f BBFF) ) ) 

1 

- (2 fAAEF + 4 fAAFF + fABEF + 2 fABFF) ( fAAEF + 2 fAAFF -(- 1 + dl ) ( fABEF + 2 fABFF) ) 
8 

1 

- (fAAEF (fABEF + 2 dl fABEF + 2 fABFF + 4 dl fABFF + 4 f BBEF + 8 f BBFF) + 
8 

2 fAAFF (fABEF + 2 dl fABEF + 2 fABFF + 4 dl fABFF + 4 f BBEF + 8 f BBFF) + 
(fABEF + 2 fABFF) ( f ABEF + 2 f ABFF - ( - 3 + 2 dl ) ( f BBEF + 2 f BBFF ) ) ) 

1 

- (dl (fABEF + 2 fABFF) + f BBEF + 2 f BBFF) (fABEF + 2 (fABFF + f BBEF + 2 f BBFF) ) 



Selection 

wAAEE = wAAEF = wAAFF = 1 ; 
wABEE = wABEF = wABFF = 1 - hs ; 
wBBEE = wBBEF = wBBFF = 1 - S ; 

W = FullSimplify [ (f AAEE Drive wAAEE + f AAEF Drive wAAEF + f AAFF Drive wAAFF 
+ f ABEE Drive wABEE + f ABEF Drive wABEF + f ABFF Drive wABFF 
+ f BBEE Drive wBBEE + f BBEF Drive wBBEF + f BBFF Drive wBBFF 

)] ; 



40 | MathematicaAppendix.nb 



fAAEE 


e i = FullSimplify [fAAEE Drive wAAEE / w] ; 










fABEE 


e i = FullSimplify [fABEE Drive wABEE / w] ; 










f BBEE 


ei = FullSimplify [fBBEE Drive wBBEE / w] ; 










fAAEF 


ei = FullSimplify [fAAEF Drive wAAEF / w] ; 










fABEF 


ei = FullSimplify [ f ABEF Drive wABEF / w] ; 










fBBEF 


ei = FullSimplify [fBBEF Drive wBBEF / w] ; 










fAAFF 


ei = FullSimplify [fAAFF Drive wAAFF / w] ; 










fABFF 


ei = FullSimplify [ f ABFF Drive wABFF / w] ; 










fBBFF 


ei = FullSimplify [ f BBFF Drive wBBFF / w] ; 










f Asel = 


FullSimplify [fAAEEgei + f AAEF Sel + f AAFF Se 


i + (fABEE S e 


l + f ABEFgei 


+ f ABFF Se 


) /2]; 


f B se i = 


FullSimplify [fBBEE S ei + f BBEF Se i + f BBFF Se 


l + (fABEE Se 


l + f ABEFgei 


+ f ABFF Se 


) /2] ; 


fE Se i = 


FullSimplify [fAAEEgei + f ABEE Sel + f BBEE Se 


l + (fAAEF S e 


l + f ABEFgei 


+ fBBEF S e 


) /2] ; 


f Fsel = 


FullSimplify [fAAFFgei + f ABFF Se i + f BBFF Se 


l + (fAAEFge 


l + f ABEFgei 


+ fBBEF Se 


) /2] ; 



Af A = FullSimplify [fA S ei - fA] 
AfB = FullSimplify [fB S ei - fB] 
Af E = FullSimplify [fEsei - fE] 
AfF = FullSimplify [fF Se i - fF] 



(*AfC=FullSimplify [fC S ei-fC] ;*) 

(*Check: do genotype freqs after selection sum to one?*) 
FullSimplify [f A Se i + fB Se i] 
FullSimplify [fE S ei + fF Se i] 



Analysis 

Invasion of sperm-acting drive modifier unlinked to the driver 

invasionOfUnlinkedModeif ier = 
FullSimplify [w AfF / f F / . randomHaps / . D -» 0 / . sum2one / . minormod] / . 
fF-»0(*this assumes no LD and random mating*) 

- (-1 + fB) fB ( (-1 + 2 fB) hs - fB s) e 

(*With some rearangement , 

invasionOfUnlikedModeifier= -e fA fB(fB s+(fA-fB) hs)) 
So long as sh«s, [ 

i.e. the driver does not have an underdominant effect on fitness] this is 
guaranteed to be positive when e is negative [i.e. the modifier dampens drive]*) 
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Fixation of sperm-acting drive modifier unlinked to the driver 

fixationOfUnlinkedModeifier = FullSimplify [ 

FullSimplify [w AfF / (1 - fF) /. randomHaps / . D -» 0 / . sum2one /. minormod] /. fF -» l] 
(*this assumes no LD and random mating*) 

(-1 + fB) fB (hs -2fBhs + fBs) e 

(*Thus the invasion and fixation conditions are equivalent*) 

Model 6'. Female drive depends on female genotype at 
an unlinked locus 

In addition two our one locus with two alleles, A (non-driving) and B (traditional driver), we have a 
second locus with loci, E and F. We are concerned with the invasion and fixation of the F allele, which 
acts to modify drive when it fertilizes a heterozygote. Drive of the B allele is dO, but changes to d1 when 
fertilized by an 'F' bearing sperm 



Setup 

Clear All ["Global' *"] 

f A = . ; f B = . ; f E = . ; f F = . ; 
f AE = f A f E + D ; 
fAF = fAfF -D; 
f BE = f B f E - D; 
f BF = f B f F + D; 

randomHaps = { f AAEE -» FullSimplify [ fAE fAE] , 

fABEE -» FullSimplify [2 fAE f BE] , fBBEE -» FullSimplify [ f BE f BE] , 

fAAEF -» FullSimplify [2 fAE fAF] , fABEF -» FullSimplify [2 (fAE f BF + fAF f BE) ] , 

fBBEF -» FullSimplify [2 f BE f BF] , fAAFF -» FullSimplify [ fAF fAF] , 

fABFF -» FullSimplify [2 fAF f BF] , fBBFF -> FullSimplify [ fBF f BF] } ; 

sum2one = { f A -» 1 - f B, fE -» 1 - f F} 

minormod = {dl-»dO + e} 
(♦assuming the sperm acting modifier additively increases drive by epsilon*) ; 

{fA -> 1 - fB, f E -> 1 - fF} 

Drive 

fAAEE Drive = FullSimplify [ 

( f AAEE * f AAEE * 1) + (f AAEE * fABEE * 1 / 2) + (f AAEE * fBBEE * 0) + ( f AAEE * fAAEF * 1 / 2) + 
(f AAEE * fABEF * 1 / 4) + (f AAEE * fBBEF * 0) + ( f AAEE * fAAFF * 0) + (f AAEE * fABFF * 0) + 
(f AAEE * fBBFF * 0) + (fABEE * f AAEE * (1 - dO) ) + (fABEE * fABEE * ( (1 - dO) / 2) ) + 
(fABEE * fBBEE * 0) + (fABEE * fAAEF * ( (1 - dO) / 2) ) + (fABEE * fABEF * ( (1 - dO) / 4) ) + 
(fABEE * fBBEF * 0) + (fABEE * fAAFF * 0) + (fABEE * fABFF * 0) + (fABEE * fBBFF * 0) + 
(fBBEE * f AAEE * 0) + (fBBEE * fABEE * 0) + (fBBEE * fBBEE * 0) + (fBBEE * fAAEF * 0) + 
(fBBEE * fABEF * 0) + (fBBEE * fBBEF * 0) + (fBBEE * fAAFF * 0) + (fBBEE * fABFF * 0) + 
(fBBEE * fBBFF * 0) + (fAAEF * f AAEE * 1 / 2) + (fAAEF * fABEE * 1 / 4) + 
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(fAAEF * fBBEE * 0) + (fAAEF * fAAEF * 1 / 4) + (fAAEF * fABEF * 1 / 8) + 
(fAAEF * fBBEF * 0) + (fAAEF * fAAFF * 0) + (fAAEF * fABFF * 0) + (fAAEF * fBBFF * 0) + 
(fABEF * fAAEE * ( (1 - dh) / 2) ) + (fABEF * fABEE * ( ( 1 - dh) / 4) ) + (fABEF * fBBEE * 0) + 
(fABEF * fAAEF * ( (1 - dh) / 4) ) + (fABEF * fABEF * ( ( 1 - dh) / 8) ) + (fABEF * fBBEF * 0) + 
(fABEF * fAAFF * 0) + (fABEF * fABFF * 0) + (fABEF * fBBFF * 0) + (fBBEF * fAAEE * 0) + 
(fBBEF * fABEE * 0) + (fBBEF * fBBEE * 0) + (fBBEF * fAAEF * 0) + (fBBEF * fABEF * 0) + 
(fBBEF * fBBEF * 0) + (fBBEF * fAAFF * 0) + (fBBEF * fABFF * 0) + (fBBEF * fBBFF * 0) + 
(fAAFF * fAAEE * 0) + (fAAFF * fABEE * 0) + (fAAFF * fBBEE * 0) + (fAAFF * fAAEF * 0) + 
(fAAFF * fABEF * 0) + (fAAFF * fBBEF * 0) + (fAAFF * fAAFF * 0) + (fAAFF * fABFF * 0) + 
(fAAFF * fBBFF * 0) + (fABFF * fAAEE * 0) + (fABFF * fABEE * 0) + (fABFF * fBBEE * 0) + 
(fABFF * fAAEF * 0) + (fABFF * fABEF * 0) + (fABFF * fBBEF * 0) + (fABFF * fAAFF * 0) + 
(fABFF * fABFF * 0) + (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 0) + (fBBFF * fABEE * 0) + 
(fBBFF * fBBEE * 0) + (fBBFF * fAAEF * 0) + (fBBFF * fABEF * 0) + (fBBFF * fBBEF * 0) + 
(fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) +0] 
fABEE Drive = FullSimplify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 1 / 2) + (fAAEE * fBBEE * 1) 
(fAAEE * fAAEF * 0) + (fAAEE * fABEF * 1 / 4) + (fAAEE * fBBEF * 1 / 2) + 

(fAAEE * fAAFF * 0) + (fAAEE * fABFF * 0) + (fAAEE * fBBFF * 0) + (fABEE * fAAEE * (dO) ) + 
(fABEE * fABEE * (1 / 2) ) + (fABEE * fBBEE * (1 - dO) ) + (fABEE * fAAEF * (dO / 2) ) + 
(fABEE * fABEF * (1 / 4) ) + (fABEE * fBBEF * ( (1 - dO) / 2) ) + (fABEE * fAAFF * 0) + 
(fABEE * fABFF * 0) + (fABEE * fBBFF * 0) + (fBBEE * fAAEE * 1) + (fBBEE * fABEE * 1 / 2) + 
(fBBEE * fBBEE * 0) + (fBBEE * fAAEF * 1 / 2) + (fBBEE * fABEF * 1 / 4) + (fBBEE * fBBEF * 0) 
(fBBEE * fAAFF * 0) + (fBBEE * fABFF * 0) + (fBBEE * fBBFF * 0) + (fAAEF * fAAEE * 0) + 
(fAAEF * fABEE * 1 / 4) + (fAAEF * fBBEE * 1 / 2) + (fAAEF * fAAEF * 0) + 

(fAAEF * fABEF * 1 / 8) + (fAAEF * fBBEF * 1 / 4) + (fAAEF * fAAFF * 0) + (fAAEF * fABFF * 0) 
(fAAEF * fBBFF * 0) + (fABEF * fAAEE * (dh / 2) ) + (fABEF * fABEE * (1 / 4) ) + 
(fABEF * fBBEE * ( (1 - dh) / 2) ) + (fABEF * fAAEF * (dh / 4) ) + (fABEF * fABEF * (1 / 8) ) + 
(fABEF * fBBEF * ( (1 - dh) / 4) ) + (fABEF * fAAFF * 0) + (fABEF * fABFF * 0) + 
(fABEF * fBBFF * 0) + (fBBEF * fAAEE * 1 / 2) + (fBBEF * fABEE * 1 / 4) + 
(fBBEF * fBBEE * 0) + (fBBEF * fAAEF * 1 / 4) + (fBBEF * fABEF * 1 / 8) + 
(fBBEF * fBBEF * 0) + (fBBEF * fAAFF * 0) + (fBBEF * fABFF * 0) + (fBBEF * fBBFF * 0) + 
(fAAFF * fAAEE * 0) + (fAAFF * fABEE * 0) + (fAAFF * fBBEE * 0) + (fAAFF * fAAEF * 0) + 
(fAAFF * fABEF * 0) + (fAAFF * fBBEF * 0) + (fAAFF * fAAFF * 0) + (fAAFF * fABFF * 0) + 
(fAAFF * fBBFF * 0) + (fABFF * fAAEE * 0) + (fABFF * fABEE * 0) + (fABFF * fBBEE * 0) + 
(fABFF * fAAEF * 0) + (fABFF * fABEF * 0) + (fABFF * fBBEF * 0) + (fABFF * fAAFF * 0) + 
(fABFF * fABFF * 0) + (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 0) + (fBBFF * fABEE * 0) + 
(fBBFF * fBBEE * 0) + (fBBFF * fAAEF * 0) + (fBBFF * fABEF * 0) + (fBBFF * fBBEF * 0) + 
(fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) +0] 
fBBEE Drive = FullSimplify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
(fAAEE * fAAEF * 0) + (fAAEE * fABEF * 0) + (fAAEE * fBBEF * 0) + (fAAEE * fAAFF * 0) + 
(fAAEE * fABFF * 0) + (fAAEE * fBBFF * 0) + (fABEE * fAAEE * 0) + (fABEE * fABEE * (dO / 2) ) 
(fABEE * fBBEE * (dO) ) + (fABEE * fAAEF * 0) + (fABEE * fABEF * (dO / 4) ) + 
(fABEE * fBBEF * (dO / 2) ) + (fABEE * fAAFF * 0) + (fABEE * fABFF * 0) + (fABEE * fBBFF * 0) 
(fBBEE * fAAEE * 0) + (fBBEE * fABEE * 1 / 2) + (fBBEE * fBBEE * 1) + (fBBEE * fAAEF * 0) + 
(fBBEE * fABEF * 1 / 4) + (fBBEE * fBBEF * 1 / 2) + (fBBEE * fAAFF * 0) + (fBBEE * fABFF * 0) 
(fBBEE * fBBFF * 0) + (fAAEF * fAAEE * 0) + (fAAEF * fABEE * 0) + (fAAEF * fBBEE * 0) + 
(fAAEF * fAAEF * 0) + (fAAEF * fABEF * 0) + (fAAEF * fBBEF * 0) + (fAAEF * fAAFF * 0) + 
(fAAEF * fABFF * 0) + (fAAEF * fBBFF * 0) + (fABEF * fAAEE * 0) + (fABEF * fABEE * (dh / 4) ) 
(fABEF * fBBEE * (dh / 2) ) + (fABEF * fAAEF * 0) + (fABEF * fABEF * (dh / 8) ) + 
(fABEF * fBBEF * (dh / 4) ) + (fABEF * fAAFF * 0) + (fABEF * fABFF * 0) + (fABEF * fBBFF * 0) 
(fBBEF * fAAEE * 0) + (fBBEF * fABEE * 1 / 4) + (fBBEF * fBBEE * 1 / 2) + (fBBEF * fAAEF * 0) 
(fBBEF * fABEF * 1 / 8) + (fBBEF * fBBEF * 1 / 4) + (fBBEF * fAAFF * 0) + (fBBEF * fABFF * 0) 
(fBBEF * fBBFF * 0) + (fAAFF * fAAEE * 0) + (fAAFF * fABEE * 0) + (fAAFF * fBBEE * 0) + 
(fAAFF * fAAEF * 0) + (fAAFF * fABEF * 0) + (fAAFF * fBBEF * 0) + (fAAFF * fAAFF * 0) + 
(fAAFF * fABFF * 0) + (fAAFF * fBBFF * 0) + (fABFF * fAAEE * 0) + (fABFF * fABEE * 0) + 
(fABFF * fBBEE * 0) + (fABFF * fAAEF * 0) + (fABFF * fABEF * 0) + (fABFF * fBBEF * 0) + 
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(fABFF 
(fBBFF 
(fBBFF 
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(fBBEE 
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(fAAEF 
(fAAEF 
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(fABEF 
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fBBEF Drive = 
(fAAEE • 
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(fABEE 



FullSimp: 

fAAEF 

fABFF 



fBBFF * 0 

fAAEF 

fABFF 

fBBEE 

fBBEF 

fBBFF 

fBBEE 

fBBEF 

fBBFF 

fAAEF 

fABFF 

fBBEE 

fBBEF 

fAAEE 

fAAEF 

fAAFF 

fABEE 

fBBEF 

FullSimp. 

fAAEF 

fABFF 

fBBEE 

fBBEF 

fBBFF 

fAAEF * 1 

fABFF * 1 

fBBEE * 1 

fBBEF * 1 

fBBFF * 1 

fBBEE * ( 

fBBEF * ( 

fBBFF * ( 

fBBEE * 0 

fBBEF * 0 

fBBFF * 0 

fAAEF * 0 

fAAFF * 0 

fABEE * ( 

fABEF * ( 

fABFF * 0 

fBBEE * 0 

fBBEF * 0 

FullSimp: 

fAAEF 

fABFF 

fBBEE 



+ (fABFF * fABFF * 0) + (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 0) + 
+ (fBBFF * fBBEE * 0) + (fBBFF * fAAEF * 0) + (fBBFF * fABEF * 0) + 
+ (fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) +0] 
.ify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
2) + (fAAEE * fABEF * 1 / 4) + (fAAEE * fBBEF * 0) + (fAAEE * fAAFF * 1) 
2) + (fAAEE * fBBFF * 0) + (fABEE * fAAEE * 0) + (fABEE * fABEE * 0) + 
+ (fABEE * fAAEF * ((1-dO) / 2) ) + (fABEE * fABEF * ((1-dO) / 4) ) + 
+ (fABEE * fAAFF * (1 - dO) ) + (fABEE * fABFF * ( (1 - dO) / 2) ) + 
+ (fBBEE * fAAEE * 0) + (fBBEE * fABEE * 0) + (fBBEE * fBBEE * 0) + 
+ (fBBEE * fABEF * 0) + (fBBEE * fBBEF * 0) + (fBBEE * fAAFF * 0) + 
+ (fBBEE * fBBFF * 0) + (fAAEF * fAAEE * 1 / 2) + (fAAEF * fABEE * 1 / 4) 
+ (fAAEF * fAAEF * 1 / 2) + (fAAEF * fABEF * 1 / 4) + 
+ (fAAEF * fAAFF * 1 / 2) + (fAAEF * fABFF * 1 / 4) + 

+ (fABEF * fAAEE * ( ( 1 - dh) / 2) ) + (fABEF * fABEE * ( (1 - dh) / 4) ) + 
+ (fABEF * fAAEF * ( ( 1 - dh) / 2) ) + (fABEF * fABEF * ( (1 - dh) / 4) ) + 
+ (fABEF * fAAFF * ((1-dh) / 2) ) + (fABEF * fABFF * ((1-dh) / 4) ) + 
+ (fBBEF * fAAEE * 0) + (fBBEF * fABEE * 0) + (fBBEF * fBBEE * 0) + 
+ (fBBEF * fABEF * 0) + (fBBEF * fBBEF * 0) + (fBBEF * fAAFF * 0) + 
+ (fBBEF * fBBFF * 0) + (fAAFF * fAAEE * 1) + (fAAFF * fABEE * 1 / 2) + 
+ (fAAFF * fAAEF * 1 / 2) + (fAAFF * fABEF * 1 / 4) + 
+ (fAAFF * fAAFF * 0) + (fAAFF * fABFF * 0) + (fAAFF * fBBFF * 0) + 
-dl) ) + (fABFF * fABEE * ( (1 - dl) / 2) ) + (fABFF * fBBEE * 0) + 
1 - dl) / 2) ) + (fABFF * fABEF * ( (1 - dl) / 4) ) + (fABFF * fBBEF * 0) + 
+ (fABFF * fABFF * 0) + (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 0) + 
+ (fBBFF * fBBEE * 0) + (fBBFF * fAAEF * 0) + (fBBFF * fABEF * 0) + 
+ (fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) +0] 

.ify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
+ (fAAEE * fABEF * 1 / 4) + (fAAEE * fBBEF * 1 / 2) + (fAAEE * fAAFF * 0) 

/ 2) + (fAAEE * fBBFF * 1) + (fABEE * fAAEE * 0) + (fABEE * fABEE * 0) + 
+ (fABEE * fAAEF * (dO / 2) ) + (fABEE * fABEF * (1 / 4) ) + 
1 - dO) / 2) ) + (fABEE * fAAFF * (dO) ) + (fABEE * fABFF * (1 / 2) ) + 
- dO) ) + (fBBEE * fAAEE * 0) + (fBBEE * fABEE * 0) + (fBBEE * fBBEE * 0) 
2) + (fBBEE * fABEF * 1 / 4) + (fBBEE * fBBEF * 0) + (fBBEE * fAAFF * 1) 
2) + (fBBEE * fBBFF * 0) + (fAAEF * fAAEE * 0) + (fAAEF * fABEE * 1 / 4) 
2) + (fAAEF * fAAEF * 0) + (fAAEF * fABEF * 1 / 4) + 
2) + (fAAEF * fAAFF * 0) + (fAAEF * fABFF * 1 / 4) + 
2) + (fABEF * fAAEE * (dh / 2) ) + (fABEF * fABEE * (1 / 4) ) + 

;i - dh) / 2) ) + (fABEF * fAAEF * (dh / 2) ) + (fABEF * fABEF * (1 / 4) ) + 

;i - dh) / 2) ) + (fABEF * fAAFF * (dh / 2) ) + (fABEF * fABFF * (1 / 4) ) + 
1 - dh) / 2) ) + (fBBEF * fAAEE * 1 / 2) + (fBBEF * fABEE * 1 / 4) + 
+ (fBBEF * fAAEF * 1 / 2) + (fBBEF * fABEF * 1 / 4) + 
+ (fBBEF * fAAFF * 1 / 2) + (fBBEF * fABFF * 1 / 4) + 

+ (fAAFF * fAAEE * 0) + (fAAFF * fABEE * 1 / 2) + (fAAFF * fBBEE * 1) + 
+ (fAAFF * fABEF * 1 / 4) + (fAAFF * fBBEF * 1 / 2) + 

+ (fAAFF * fABFF * 0) + (fAAFF * fBBFF * 0) + (fABFF * fAAEE * (dl) ) + 

/ 2) ) + (fABFF * fBBEE * (1 - dl) ) + (fABFF * fAAEF * (dl / 2) ) + 

/ 4) ) + (fABFF * fBBEF * ( (1 - dl) / 2) ) + (fABFF * fAAFF * 0) + 

+ (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 1) + (fBBFF * fABEE * 1 / 2) + 

+ (fBBFF * fAAEF * 1 / 2) + (fBBFF * fABEF * 1 / 4) + 

+ (fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) +0] 
.ify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
+ (fAAEE * fABEF * 0) + (fAAEE * fBBEF * 0) + (fAAEE * fAAFF * 0) + 
+ (fAAEE * fBBFF * 0) + (fABEE * fAAEE * 0) + (fABEE * fABEE * 0) + 
(fABEE * fAAEF * 0) + (fABEE * fABEF * (dO / 4) ) + 



fBBEF * (dO / 2) ) + (fABEE * fAAFF * 0) + (fABEE * fABFF * (dO / 2) ) 4 
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(fABEE * fBBFF * dO) + ( f BBEE * fAAEE * 0) + ( f BBEE * fABEE * 0) + (fBBEE * fBBEE * 0) + 
(fBBEE * fAAEF * 0) + (fBBEE * fABEF * 1 / 4) + (fBBEE * fBBEF * 1 / 2) + (fBBEE * fAAFF * 0) 
(fBBEE * fABFF * 1 / 2) + (fBBEE * fBBFF * 1) + (fAAEF * fAAEE * 0) + (fAAEF * fABEE * 0) + 
(fAAEF * fBBEE * 0) + (fAAEF * fAAEF * 0) + (fAAEF * fABEF * 0) + (fAAEF * fBBEF * 0) + 
(fAAEF * fAAFF * 0) + (fAAEF * fABFF * 0) + (fAAEF * fBBFF * 0) + (fABEF * fAAEE * 0) + 
(fABEF * fABEE * (dh / 4) ) + (fABEF * fBBEE * (dh / 2) ) + (fABEF * fAAEF * 0) + 
(fABEF * fABEF * (dh / 4) ) + (fABEF * fBBEF * (dh / 2) ) + (fABEF * fAAFF * 0) + 
(fABEF * fABFF * (dh / 4) ) + (fABEF * fBBFF * (dh / 2) ) + (fBBEF * fAAEE * 0) + 
(fBBEF * fABEE * 1 / 4) + (fBBEF * fBBEE * 1 / 2) + (fBBEF * fAAEF * 0) + 
(fBBEF * fABEF * 1 / 4) + (fBBEF * fBBEF * 1 / 2) + (fBBEF * fAAFF * 0) + 
(fBBEF * fABFF * 1 / 4) + (fBBEF * fBBFF * 1 / 2) + (fAAFF * fAAEE * 0) + 
(fAAFF * fABEE * 0) + (fAAFF * fBBEE * 0) + (fAAFF * fAAEF * 0) + (fAAFF * fABEF * 0) + 
(fAAFF * fBBEF * 0) + (fAAFF * fAAFF * 0) + (fAAFF * fABFF * 0) + (fAAFF * fBBFF * 0) + 
(fABFF * fAAEE * 0) + (fABFF * fABEE * (dl / 2) ) + (fABFF * fBBEE * (dl) ) + 
(fABFF * fAAEF * 0) + (fABFF * fABEF * (dl / 4) ) + (fABFF * fBBEF * (dl / 2) ) + 
(fABFF * fAAFF * 0) + (fABFF * fABFF * 0) + (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 0) + 
(fBBFF * fABEE * 1 / 2) + (fBBFF * fBBEE * 1) + (fBBFF * fAAEF * 0) + (fBBFF * fABEF * 1 / 4) 
(fBBFF * fBBEF * 1 / 2) + (fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) + 0 
fAAFF Drive = FullSimplify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
(fAAEE * fAAEF * 0) + (fAAEE * fABEF * 0) + (fAAEE * fBBEF * 0) + (fAAEE * fAAFF * 0) + 
(fAAEE * fABFF * 0) + (fAAEE * fBBFF * 0) + (fABEE * fAAEE * 0) + (fABEE * fABEE * 0) + 
(fABEE * fBBEE * 0) + (fABEE * fAAEF * 0) + (fABEE * fABEF * 0) + (fABEE * fBBEF * 0) + 
(fABEE * fAAFF * 0) + (fABEE * fABFF * 0) + (fABEE * fBBFF * 0) + (fBBEE * fAAEE * 0) + 
(fBBEE * fABEE * 0) + (fBBEE * fBBEE * 0) + (fBBEE * fAAEF * 0) + (fBBEE * fABEF * 0) + 
(fBBEE * fBBEF * 0) + (fBBEE * fAAFF * 0) + (fBBEE * fABFF * 0) + (fBBEE * fBBFF * 0) + 
(fAAEF * fAAEE * 0) + (fAAEF * fABEE * 0) + (fAAEF * fBBEE * 0) + (fAAEF * fAAEF * 1 / 4) + 
(fAAEF * fABEF * 1 / 8) + (fAAEF * fBBEF * 0) + (fAAEF * fAAFF * 1 / 2) + 

(fAAEF * fABFF * 1 / 4) + (fAAEF * fBBFF * 0) + (fABEF * fAAEE * 0) + (fABEF * fABEE * 0) + 
(fABEF * fBBEE * 0) + (fABEF * fAAEF * ( ( 1 - dh) / 4) ) + (fABEF * fABEF * ( (1 - dh) / 8) ) + 
(fABEF * fBBEF * 0) + (fABEF * fAAFF * ( ( 1 - dh) / 2) ) + (fABEF * fABFF * ( (1 - dh) / 4) ) + 
(fABEF * fBBFF * 0) + (fBBEF * fAAEE * 0) + (fBBEF * fABEE * 0) + (fBBEF * fBBEE * 0) + 
(fBBEF * fAAEF * 0) + (fBBEF * fABEF * 0) + (fBBEF * fBBEF * 0) + (fBBEF * fAAFF * 0) + 
(fBBEF * fABFF * 0) + (fBBEF * fBBFF * 0) + (fAAFF * fAAEE * 0) + (fAAFF * fABEE * 0) + 
(fAAFF * fBBEE * 0) + (fAAFF * fAAEF * 1 / 2) + (fAAFF * fABEF * 1 / 4) + (fAAFF * fBBEF * 0) 
(fAAFF * fAAFF * 1) + (fAAFF * fABFF * 1 / 2) + (fAAFF * fBBFF * 0) + (fABFF * fAAEE * 0) + 
(fABFF * fABEE * 0) + (fABFF * fBBEE * 0) + (fABFF * fAAEF * ( (1 - dl) / 2) ) + 
(fABFF * fABEF * ( (1 - dl) / 4) ) + (fABFF * fBBEF * 0) + (fABFF * fAAFF * (1 - dl) ) + 
(fABFF * fABFF * ( (1 - dl) / 2) ) + (fABFF * fBBFF * 0) + (fBBFF * fAAEE * 0) + 
(fBBFF * fABEE * 0) + (fBBFF * fBBEE * 0) + (fBBFF * fAAEF * 0) + (fBBFF * fABEF * 0) + 
(fBBFF * fBBEF * 0) + (fBBFF * fAAFF * 0) + (fBBFF * fABFF * 0) + (fBBFF * fBBFF * 0) +0] 
fABFF Drive = FullSimplify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
(fAAEE * fAAEF * 0) + (fAAEE * fABEF * 0) + (fAAEE * fBBEF * 0) + (fAAEE * fAAFF * 0) + 
(fAAEE * fABFF * 0) + (fAAEE * fBBFF * 0) + (fABEE * fAAEE * 0) + (fABEE * fABEE * 0) + 
(fABEE * fBBEE * 0) + (fABEE * fAAEF * 0) + (fABEE * fABEF * 0) + (fABEE * fBBEF * 0) + 
(fABEE * fAAFF * 0) + (fABEE * fABFF * 0) + (fABEE * fBBFF * 0) + (fBBEE * fAAEE * 0) + 
(fBBEE * fABEE * 0) + (fBBEE * fBBEE * 0) + (fBBEE * fAAEF * 0) + (fBBEE * fABEF * 0) + 
(fBBEE * fBBEF * 0) + (fBBEE * fAAFF * 0) + (fBBEE * fABFF * 0) + (fBBEE * fBBFF * 0) + 
(fAAEF * fAAEE * 0) + (fAAEF * fABEE * 0) + (fAAEF * fBBEE * 0) + (fAAEF * fAAEF * 0) + 
(fAAEF * fABEF * 1 / 8) + (fAAEF * fBBEF * 1 / 4) + (fAAEF * fAAFF * 0) + 

(fAAEF * fABFF * 1 / 4) + (fAAEF * fBBFF * 1 / 2) + (fABEF * fAAEE * 0) + (fABEF * fABEE * 0) 

(fABEF * fBBEE * 0) + (fABEF * fAAEF * (dh / 4) ) + (fABEF * fABEF * (1 / 8) ) + 

(fABEF * fBBEF * ( (1 - dh) / 4) ) + (fABEF * fAAFF * (dh / 2) ) + (fABEF * fABFF * (1 / 4) ) + 

(fABEF * fBBFF * ( (1 - dh) / 2) ) + (fBBEF * fAAEE * 0) + (fBBEF * fABEE * 0) + 

(fBBEF * fBBEE * 0) + (fBBEF * fAAEF * 1 / 4) + (fBBEF * fABEF * 1 / 8) + (fBBEF * fBBEF * 0) 

(fBBEF * fAAFF * 1 / 2) + (fBBEF * fABFF * 1 / 4) + (fBBEF * fBBFF * 0) + (fAAFF * fAAEE * 0) 
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(fAAFF * f ABEE * 0) + (fAAFF * fBBEE * 0) + (fAAFF * fAAEF * 0) + (fAAFF * fABEF * 1 / 4) + 
(fAAFF * fBBEF * 1 / 2) + (fAAFF * fAAFF * 0) + (fAAFF * fABFF * 1 / 2) + (fAAFF * fBBFF * 1) 
(fABFF * fAAEE * 0) + (fABFF * fABEE * 0) + (fABFF * fBBEE * 0) + (fABFF * fAAEF * (dl / 2) ) 
(fABFF * fABEF * (1 / 4) ) + (fABFF * fBBEF * ( (1 - dl) / 2) ) + (fABFF * fAAFF * (dl) ) + 
(fABFF * fABFF * (1 / 2) ) + (fABFF * fBBFF * (1 - dl) ) + (fBBFF * fAAEE * 0) + 
(fBBFF * fABEE * 0) + (fBBFF * fBBEE * 0) + (fBBFF * fAAEF * 1 / 2) + (fBBFF * fABEF * 1 / 4) 
(fBBFF * fBBEF * 0) + (fBBFF * fAAFF * 1) + (fBBFF * fABFF * 1 / 2) + (fBBFF * fBBFF * 0) +0] 
fBBFF Drive = FullSimplify [ (fAAEE * fAAEE * 0) + (fAAEE * fABEE * 0) + (fAAEE * fBBEE * 0) + 
(fAAEE * fAAEF * 
(fAAEE * fABFF * 
(fABEE * fBBEE * 
(fABEE * fAAFF * 
(fBBEE * fABEE * 
(fBBEE * fBBEF * 
(fAAEF * fAAEE * 
(fAAEF * fABEF * 
(fAAEF * fBBFF * 
(fABEF * fAAEF * 
(fABEF * fAAFF * 
(fBBEF * fAAEE * 

(fBBEF * fABEF * 1 / 8) + (fBBEF * fBBEF * 1 / 4) + (fBBEF * fAAFF * 0) + 

(fBBEF * fABFF * 1 / 4) + (fBBEF * fBBFF * 1 / 2) + (fAAFF * fAAEE * 0) + (fAAFF * fABEE * 0) 
(fAAFF * fBBEE * 0) + (fAAFF * fAAEF * 0) + (fAAFF * fABEF * 0) + (fAAFF * fBBEF * 0) + 
(fAAFF * fAAFF * 0) + (fAAFF * fABFF * 0) + (fAAFF * fBBFF * 0) + (fABFF * fAAEE * 0) + 
(fABFF * fABEE * 0) + (fABFF * fBBEE * 0) + (fABFF * fAAEF * 0) + (fABFF * fABEF * (dl / 4) ) 
(fABFF * fBBEF * (dl / 2) ) + (fABFF * fAAFF * 0) + (fABFF * fABFF * (dl / 2) ) + 
(fABFF * fBBFF * dl) + (fBBFF * fAAEE * 0) + (fBBFF * fABEE * 0) + (fBBFF * fBBEE * 0) + 
(fBBFF * fAAEF * 0) + (fBBFF * fABEF * 1 / 4) + (fBBFF * fBBEF * 1 / 2) + 
(fBBFF * fAAFF * 0) + (fBBFF * fABFF * 1 / 2) + (fBBFF * fBBFF * 1) +0] 

1 

- (4 fAAEE + 2 (fAAEF + fABEE) + fABEF) 
8 

(2 fAAEE + fAAEF + 2 fABEE - 2 dO fABEE + fABEF - dh fABEF) 



0) 




(fAAEE * 


fABEF 


* 0) + (fAAEE * fBBEF 


*0) 




(fAAEE 


* fAAFF 


*0) 


0) 




(fAAEE * 


fBBFF 


* 0) + (fABEE * fAAEE 


*0) 




( fABEE 


* fABEE 


*0) 


0) 




( fABEE * 


fAAEF 


* 0) + (fABEE * fABEF 


*0) 




( fABEE 


* fBBEF 


*0) 


0) 




( fABEE * 


fABFF 


* 0) + (fABEE * fBBFF 


*0) 




( fBBEE 


* fAAEE 


*0) 


0) 




( fBBEE * 


fBBEE 


* 0) + (fBBEE * fAAEF 


*0) 




( fBBEE 


* fABEF 


*0) 


0) 




( fBBEE * 


fAAFF 


* 0) + (fBBEE * fABFF 


*0) 




( fBBEE 


* fBBFF 


*0) 


0) 




(fAAEF * 


fABEE 


* 0) + (fAAEF * fBBEE 


*0) 




(fAAEF 


* fAAEF 


*0) 


0) 




(fAAEF * 


fBBEF 


* 0) + (fAAEF * fAAFF 


*0) 




(fAAEF 


* fABFF 


*0) 


0) 




( fABEF * 


fAAEE 


* 0) + (fABEF * fABEE 


*0) 




( fABEF 


* fBBEE 


*0) 


0) 




( fABEF * 


fABEF 


* (dh / 8) ) + (fABEF * 


fBBEF * (dh / 


4)) + 




0) 




( fABEF * 


fABFF 


* (dh/4)) + (fABEF* 


fBBFF * (dh / 


2)) + 




0) 




( fBBEF * 


fABEE 


* 0) + (fBBEF * fBBEE 


*0) 




( fBBEF 


* fAAEF 


*0) 



- ( 4 fABEE 2 + 4 fABEE fABEF + fABEF 2 + 12 fABEE fBBEE - 8 dO fABEE fBBEE + 

8 V 

6 fABEF fBBEE - 4 dh fABEF fBBEE + 6 fABEE fBBEF - 4 dO fABEE fBBEF + 3 fABEF fBBEF - 

2 dh fABEF fBBEF + fAAEF (2 fABEE + 4 dO fABEE + fABEF + 2 dh fABEF + 8 fBBEE + 4 fBBEF) + 

2 fAAEE ( (2 + 4 dO) fABEE + fABEF + 2 dh fABEF + 8 fBBEE + 4 fBBEF) ) 

1 

- (2 dO fABEE + dh fABEF + 2 fBBEE + fBBEF) (2 fABEE + fABEF + 4 fBBEE + 2 fBBEF) 
8 

- ( 2 fAAEF 2 + 6 fAAFF fABEE - 4 dO fAAFF fABEE + 3 fAAFF fABEF - 

4 X 

2 dh fAAFF fABEF + 2 fABEE fABEF - dO fABEE fABEF - dh fABEE fABEF + 

fABEF 2 - dh fABEF 2 - (2 (-2 + dO + dl) fABEE + (-2 + dl + dh) fABEF) fABFF + 

fAAEE (4 fAAEF + 8 fAAFF + 3 fABEF - 2 dh fABEF + 6 fABFF - 4 dl fABFF) + 

fAAEF (4 fAAFF + 3 fABEE - 2 dO fABEE + 3 fABEF - 2 dh fABEF + 3 fABFF - 2 dl fABFF) ) 
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- (2 f ABEE fABEF + fABEF 2 + 4 fABEE fABFF + 2 fABEF fABFF + 3 fABEF fBBEE - 
4 

2 dh fABEF fBBEE + 6 fABFF fBBEE - 4 dl fABFF fBBEE + 3 fABEE f BBEF - 
2 dO fABEE f BBEF + 3 fABEF f BBEF - 2 dh fABEF f BBEF + 3 fABFF f BBEF - 
2 dl fABFF f BBEF + fAAFF ( (2 + 4 dO) fABEE + fABEF + 2 dh fABEF + 8 fBBEE + 4 f BBEF) + 
6 fABEE f BBFF - 4 dO fABEE f BBFF + 3 fABEF f BBFF - 2 dh fABEF f BBFF + 
fAAEE ( fABEF + 2 dh fABEF + 2 fABFF + 4 dl fABFF + 4 f BBEF + 8 f BBFF ) + fAAEF 
( fABEE + 2 dO fABEE + fABEF + 2 dh fABEF + fABFF + 2 dl fABFF + 4 ( fBBEE + f BBEF + f BBFF ) ) ) 

1 

- (2 dl fABEE fABFF + dl fABEF fABFF + fABEF fBBEE + 2 fABFF fBBEE + 
4 

4 dl fABFF fBBEE + fABEE f BBEF + fABEF f BBEF + fABFF f BBEF + 2 dl fABFF f BBEF + 
4 fBBEE f BBEF + 2 f BBEF 2 + (2 fABEE + fABEF + 8 fBBEE + 4 f BBEF ) f BBFF + 
dh fABEF (fABEE + fABEF + fABFF + 2 (fBBEE + f BBEF + f BBFF) ) + 
dO fABEE ( fABEF + 2 ( fABFF + f BBEF + 2 f BBFF) ) ) 

1 

- (2 fAAEF + 4 fAAFF + fABEF + 2 fABFF) 
8 

(fAAEF + 2 fAAFF + fABEF - dh fABEF + 2 fABFF - 2 dl fABFF) 

1 

- ( fABEF 2 + 4 fABEF fABFF + 4 fABFF 2 + 3 fABEF f BBEF - 2 dh fABEF f BBEF + 6 fABFF f BBEF - 

8 V 

4 dl fABFF f BBEF + 6 fABEF f BBFF - 4 dh fABEF f BBFF + 12 fABFF f BBFF - 

8 dl fABFF f BBFF + fAAEF (fABEF + 2 dh fABEF + 2 fABFF + 4 dl fABFF + 4 f BBEF + 8 f BBFF ) + 
2 fAAFF (fABEF + 2 dh fABEF + 2 fABFF + 4 dl fABFF + 4 f BBEF + 8 f BBFF) ) 

1 

- (dh fABEF + 2 dl fABFF + f BBEF + 2 f BBFF) (fABEF + 2 (fABFF + f BBEF + 2 f BBFF) ) 



Selection 

wAAEE = wAAEF = wAAFF = 1 ; 
wABEE = wABEF = wABFF = 1 - hs ; 
wBBEE = wBBEF = wBBFF = 1 - S ; 

W = FullSimplify [ (f AAEE Drive wAAEE + f AAEF Drive wAAEF + f AAFF Drive wAAFF 
+ f ABEE Drive wABEE + f ABEF Drive wABEF + f ABFF Drive wABFF 
+ f BBEE Drive wBBEE + f BBEF Drive wBBEF + f BBFF Drive wBBFF 

)] ; 
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fAAEE 


e i = FullSimplify [fAAEE Drive wAAEE / w] ; 










fABEE 


e i = FullSimplify [fABEE Drive wABEE / w] ; 










fBBEE 


ei = FullSimplify [fBBEE Drive wBBEE / w] ; 










fAAEF 


ei = FullSimplify [fAAEF Drive wAAEF / w] ; 










fABEF 


ei = FullSimplify [ f ABEF Drive wABEF / w] ; 










fBBEF 


ei = FullSimplify [fBBEF Drive wBBEF / w] ; 










fAAFF 


ei = FullSimplify [fAAFF Drive wAAFF / w] ; 










fABFF 


ei = FullSimplify [ f ABFF Drive wABFF / w] ; 










fBBFF 


ei = FullSimplify [ f BBFF Drive wBBFF / w] ; 










f Asel = 


FullSimplify [fAAEEgei + f AAEF Sel + f AAFF Se 


i + (fABEE S e 


l + f ABEFgei 


+ f ABFF Se 


) /2]; 


f B se i = 


FullSimplify [fBBEE S ei + f BBEF Se i + f BBFF Se 


l + (fABEE Se 


l + f ABEFgei 


+ f ABFF Se 


) /2] ; 


fE Se i = 


FullSimplify [fAAEEgei + f ABEE Sel + f BBEE Se 


l + (fAAEF S e 


l + f ABEFgei 


+ fBBEF S e 


) /2] ; 


f Fsel = 


FullSimplify [fAAFFgei + f ABFF Se i + f BBFF Se 


l + (fAAEFge 


l + f ABEFgei 


+ fBBEF Se 


) /2] ; 



Af A = FullSimplify [fA S ei - fA] 
AfB = FullSimplify [fB S ei - fB] 
Af E = FullSimplify [fEsei - fE] 
AfF = FullSimplify [fF Se i - fF] 



(*AfC=FullSimplify [fC S ei-fC] ;*) 

(*Check: do genotype freqs after selection sum to one?*) 
FullSimplify [f A Se i + fB Se i] 
FullSimplify [fE S ei + fF Se i] 



Analysis 

Invasion of egg-acting drive modifier unlinked to the driver 

invasionOfUnlinkedModifier = 
FullSimplify [w AfF / f F / . randomHaps / . D -» 0 / . dh -» dl / . sum2one / . minormod] / . 
fF-»0(*this assumes no LD and random mating*) 

- (-1 + fB) fB ( (-1 + 2 fB) hs - fB s) e 

(*With some rearangement , 

invasionOfUnlikedModeifier= -e fA fB(fB s+(fA-fB) hs)) 
So long as sh«s, [ 

i.e. the driver does not have an underdominant effect on fitness] this is 
guaranteed to be positive when e is negative [i.e. the modifier dampens drive]*) 
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Fixation of eggs-acting drive modifier unlinked to the driver 

fixationOfUnlinkedModifier = FullSimplif y [FullSimplify [ 

W AfF / (1 - fF) ^2 / . randomHaps / .D-»0/.dh-»dl/. sum2one / . minormod / . 
hs -» 0] ] /. fF->l(*this assumes no LD and random mating*) 

(-1 + fB) fB 2 s e 

(*Thus the invasion and fixation conditions are equivalent*) 



